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Abstract 

The lowest energy levels of positronium are studied in the non- 
relativistic approximation within the framework of Relativistic Schro- 
dinger Theory (RST). Since it is very difficult to find the exact solu- 
tions of the RST field equations (even in the non-relativistic limit), an 
approximation scheme is set up on the basis of the hydrogen-like wave 
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functions (i.e. polynomial times exponential). For any approximation 
order N (N = 0,1,2,3,...) there arises a spectrum of approximate 
RST solutions with the associated energies, quite similarly to the con- 
ventional treatment of positronium in the standard quantum theory 
(Appendix). For the lowest approximation order (N = 0) the RST 
prediction for the groundstate energy exactly agrees with the con- 
ventional prediction of the standard theory. However for the higher 
approximation orders (N = 1, 2, 3), the corresponding RST prediction 
differs from the conventional result by (roughly) 0, 9 [eV\ which con- 
firms the previous estimate of the error being due to the use of the 
spherically symmetric approximation. The excited states require the 
application of higher-order approximations (N >> 3) and are there- 
fore not adequately described by the present orders (N < 3). 



PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 
03.65.Ge - Solutions of Wave Equations: Bound States; 
03.65.Sq - Semiclassical Theories and Applications; 03.75.B 

- Matter Waves 
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I Introduction and Survey of Results 



Apparently the general belief says that the standard quantum theory is 
unique in the sense that no other theory can quantitatively predict the phe- 
nomena of the microscopic world more accurately than just that conventional 
theory. But this does logically not exclude the possibility, that other forms 
of quantum theory may perhaps exist which have the same predictive power: 
"The natural meaning of the claim that quantum theory provides for a com- 
plete scientific account of atomic phenomena is that no theoretical construc- 
tion can yield experimentally verifiable predictions about atomic phenomena 
that cannot be extracted from a quantum theoretical description" pLj. It is 
true, this claim is surely intended to refer to the standard form of quantum 
theory; however, if one could find an alternative formalism with the same 
predictive potentiality as the standard theory, this would be nevertheless of 
considerable relevance because it could shed new light on the old interpreta- 
tion problems of the standard theory. 

Perhaps the recently established Relativistic Schrodinger Theory (RST) 
does represent such an alternative candidate E]. This theory is not of 
probabilistic but rather of fluid-dynamic character and therefore resembles 
rather the well-known density functional theory |H [5] than the usual statis- 
tical approach. Naturally, such a philosophical difference must then entail 
also certain consequences for the mathematical structure of both theoretical 
approaches: the A^-particle systems are described in RST by the Whitney 
sum of the one-particle bundles, not by the tensor product of the one-particle 
Hilbert spaces as in the standard theory. However it is well-known that also 
such a fluid-dynamic approach as the density functional theory can predict 
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atomic and molecular data with the same accuracy as the standard proba- 
bilistic approach, and this should justify the endeavors for further elaborating 
those fluid-dynamic theories. 

For instance, it has recently been observed that RST, when being ap- 
plied to the non-relativistic positronium problem, predicts exactly the same 
groundstate energy as the conventional quantum mechanics [3]. Recall here 
that this conventional prediction is based upon a very simple argument: As- 
suming (for the non-relativistic situation) that the electron and positron un- 
dergo the electromagnetic interaction via the usual Coulomb potential, the 
corresponding two-particle Hamiltonian H reads 

H=pL+pt__jL_. (1.1) 

2M 2M ||ri-f2|| ^ ^ 

This Hamiltonian can easily be transformed to the sum of an external part 
(due to the center-of-mass motion) and an internal part Hi^ |6] 

An = - - (1.2) 

2m r 

which then yields the well-known conventional spectrum of the internal ex- 
citations: 

1 6,8029... , , 

AttB {rip + ly [Up + ly 

with the principal quantum number Hp being some integer (rip = 0, 1, 2, 3 . . .) 
and Ob the Bohr radius 

Indeed, the eigenvalue equation of the internal Hamiltonian fll.2p looks very 
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simple 



Hin(f) 



(m 




— ... reduced mass 
2 



T^A H ] (f) = Econv(j) 



r 



(1.5) 



and can easily be solved by the product of a polynomial times an exponential 
(cf. Appendix). 

But the problem for RST consists now just in that exact coincidence of 
its groundstate prediction [3l [7j with the conventional result fll.3l) for rip = 0! 
Namely, whereas that conventional result owns the status of absolute ex- 
actness within the standard framework (apart from relativistic and QED 
corrections [8]), its RST counterpart is only of approximative character, even 
within the non-relativistic approximation of RST. The point here is that even 
the non-relativistic approximation of RST does not lead to a spherically- 
symmetric potential (such as the conventional Coulomb potential occurring 
in the internal Hamiltonian fll.51) ): but rather the non-relativistic RST poten- 
tial must be further truncated to the spherically-symmetric form in order that 
the corresponding groundstate prediction can exactly agree with the conven- 
tional result (II. 3P for Up = 0. But even if such a truncation to the spherical 
symmetry of the RST groundstate configuration is accepted, the correspond- 
ing eigenvalue problem cannot be solved exactly; but rather one has to resort 
to some approximative procedure whose first step then yields the mentioned 
coincidence of the RST and conventional predictions, see refs. [3l [7] . 

This circumstance, i.e. the numerical coincidence of the "exact" conven- 
tional result and the approximative RST result in the non-relativistic domain, 
is just the reason why one cannot be completely satisfied with the present 
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situation. Namely, if one wishes to stick to the claim that the true nonrela- 
tivistic limit of RST should agree to a high precision with the conventional 
predictions, then the spherically symmetric RST approximation of that true 
limit should differ by a certain finite amount from both the conventional 
results and the true RST prediction. This difference, being due to the spher- 
ically symmetric approximation of the non-relativistic interaction potential, 
has been estimated as (roughly) 1 [eV], see the appendix of ref. [7j. There- 
fore one has to expect that the spherically symmetric (non-relativistic) RST 
limit of the positronium groundstate should differ from the conventional pre- 
diction fll.3|) by roughly 1 [eV]. The present paper aims at the verification 
of this energy difference and thus upholds the possibility that the true (i.e. 
anisotropic) non-relativistic RST limit coincides practically with the conven- 
tional predictions. (The common difference of both theoretical approaches 
relative to the experimental data would then refer exclusively to the rela- 
tivistic corrections). 

In this sense, the main result of the present paper consists in the elabora- 
tion of the expected energy difference of roughly 1 [eV]. This goal is attained 
by setting up a rigorous approximation formalism so that the incidental nu- 
merical coincidence of the RST and conventional predictions do appear only 
in the zero-order approximation. In the first-order approximation (N = 1), 
there emerges an energy difference of (roughly) 0, 9 [eV] which thus verifies 
the expectations. Since this energy difference persists also in the higher- 
order approximations (N = 2, 3) it seems justified to hypothesize that the 
true non-relativistic RST limit (which includes also the anisotropic effects) 
will actually meet with the corresponding non-relativistic predictions of the 
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conventional theory. 

These results are worked out in the following arrangement: 
First, a brief sketch of the general relativistic two-particle theory is pre- 
sented in Sect. II. Here the principal RST dynamics is displayed, including 
the specialization to the stationary bound states. As a preparation of the 
subsequent treatment of the positronium groundstate, the relativistic field 
equations are then transcribed to their non-relativistic and spherically sym- 
metric approximations. 

Next, the RST principle of minimal energy [3] is presented in Sect. III. 
This variational principle is the crucial point for the subsequent calculation 
of the positronium groundstate energy since the exact solutions of the RST 
field equations are very hard to obtain, so that one has to be satisfied with 
such an approximate variational technique. But fortunately this is most 
adequately provided by that principle of minimal energy whose variational 
equations consist in the coupled matter and gauge field equations (see the 
Dirac equations flll.45al) - flll.45bl) for the matter subsystem and the Poisson 
equations flII.43aP - flII.43d|) for the gauge field subsystem). Both subsystem 
variations do need a constraint, i.e. the wave function normalization for the 
matter fields QUI. lip and the Poisson identities fllll.lQp for the gauge fields. 
The latter identities refer to the numerical equality of the pure gauge field 
energy and the interaction energy M'^'^^c^ of the matter and gauge fields 
{mass equivalent) . The non-relativistic version of the matter part of this vari- 
ational principle turns out as nothing else than the well-known Ritz principle 
of conventional quantum theory. The Poisson identities do not only ensure 
the deduction of the gauge field equations from the principle of minimal en- 
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ergy, but provide also a nice consistency check for the concrete calculation 
of the electric gauge field energy -E^"*. Namely, the latter energy essentially 
consists of three parts: the groundstate contribution, the quadratic terms 
and the quartic terms, see equation flIII.53p . Since the same subdivision does 
hold also for the electric mass equivalent M^^^c^ (IIII.6ip . the Poisson identity 
implies the separate numerical equality of all three contributions, see equa- 
tions (1III.62|) - (1III.67|) . Clearly such a set of integral identities is a rigorous 
consistency test for the correct computation of the gauge field energy with 
respect to the chosen trial functions. 

With all these preparations, the proper goal of the paper can be tackled in 
Sect. IV: Selecting as trial functions the hydrogen-like wave functions R{r) 
(see equation flIII.35p together with (lIV.ip ). one can study the action of the 
non-relativistic energy functional e!^"^ ( 1III.32I) on this set of functions, of 
course in order to determine the stationary points and especially the min- 
imally possible energy value (groundstate). Since the hydrogen-like wave 
functions are the product of an exponential times a polynomial of degree N, 
these trial functions may be classified just by the degree N. In this way 
there arises an energy spectrum for any order N (N = 0,1,2,3,...) each 
of which has a finite number of members. The lowest-energy state is of 
course the groundstate which is our main concern because its RST energy 
has been found to agree exactly with the conventional prediction (II. 3|) . In 
order to better understand this strange agreement of the RST and conven- 
tional groundstate predictions, the same variational technique is applied in 
the appendix also to the conventional positronium problem fll.5l) : and the 
comparison of both approaches yields the following results: 
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The lowest-order spectrum (N = 0) consists exclusively of the ground- 
state whose energy is found in this order to be the same for both approaches, 
cf. the RST result flIV.36|) and the conventional result flA.18|) . For the higher- 
order spectra (N > 1) it is important to observe that the constraint of wave 
function normalization forces one to look for the stationary points of the en- 
ergy function on an N-dimensional sphere S*^. This is not quite trivial from 
the topological point of view since the well-known Morse inequalities [9l [TO] 
imply certain restrictions upon the number of stationary points which a con- 
tinuously differentiable function can own on a given compact manifold. For 
instance, on an N-dimensional torus = x x . . . the minimally pos- 
sible number of stationary points is 2^ [11]. But fortunately the proposition 
of Reeb [T2| ensures for the spheres 5*^ that the minimally possible number 
is 2; and therefore the conventional approach for N = 1 Appendix) can 
produce two stationary points (i.e. minimum and maximum) of the energy 
function on the circle S^, see fig. A.l. Clearly this corresponds to the ground- 
state and first excited state as exact solutions of the internal Schrodinger 
equation fll.Sp . Furthermore, the Reeb proposition admits that the general 
conventional energy function flA.14p has exactly N+1 stationary points on S*^ 
which generate the conventional spectrum fll.Sp with < Up < N; see the 
exemplification for N = 2 in the Appendix. 

A very pleasant feature of the conventional spectrum refers to the fact, 
that all stationary points of £'conv{f^} (1A.14p on 5^ are present also on S^'^^, 
i.e. the property of stationarity is not lost by the embedding 5*^ 
This is in contrast to the situation with RST. For instance, the RST ground- 
state energy pV.36p for N = actually emerges also in the next order N = 1 
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in form of the field configuration with energy i?ia (lIV.45ap . but this corre- 
sponds in this order N = 1 to a maximum (or saddle point, resp.) and not 
to a minimum as in the lowest order N = 0. Therefore the agreement of the 
conventional result with the RST groundstate energy flIV.36p and flIV.45ap 
is not "stable" within the RST approach, as is the case within the conven- 
tional approach. The concept of "stability" refers here to the property of a 
stationary point to be preserved by the embedding S^^^. The insta- 

bility of the "exact" RST energy (1IV.36P and (lIV.45aP for the positronium 
groundstate is revealed by the next higher embedding S*^ ^ S*^. Indeed, this 
embedding does preserve three of the four stationary points of order N = 1, 
but the "exact" configuration with i?ia flIV.45aP is lost and there arise six 
new stationary points with different energies i?2e — -S2j f HV.57aP - flIV.57fp . 

However concerning the groundstate, this situation just meets with the 
expectations concerning the numerical accuracy of the present spherically 
symmetric approximation. Since it has been estimated that this kind of 
approximation generates a deviation of (roughly) 1 [eV] from the true (but 
non-relativistic) RST predictions for the positronium groundstate [7J, it is 
necessary that the present RST prediction deviates from the conventional 
result flA.lSp by roughly 1 [eV]. Now the minimal value of the RST energy 
function e!^'^ (]IV.24p is found here in the order N = 1 as Eid (]IV.46bp . and 
this survives also the passage to the next order N = 2. It is true, the case with 
the "exact" RST configuration (la) flIV.45aP demonstrates that a stationary 
point of a lower order may survive to some higher order but nevertheless 
vanishes during the passage to the subsequent orders, but the relatively small 
change of the RST groundstate energy for N = (-^ Eq = —6, 8029 . . . [eV]), 
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cf. firOej) . to its value for N = 1 and N = 2 = -7, 675 . . . [eV]) 

suggests to hypothesize that the "true" (i.e. N oo) RST prediction (in 
the spherically-symmetric approximation) actually remains sufficiently close 
to the exact conventional value fllV.36P : i.e. more precisely: the isotropic 
RST prediction remains within the limits which are put by the anisotropic 
estimate [Zj of 1 [eV], so that the "true" RST prediction (i.e. with observation 
of the anisotropy) practically can agree with the conventional result! Clearly 
in order to verify this rigorously, one would have to go explicitly through 
the higher orders up to N — > oo. Here one could then also study the more 
mathematical questions of how many and which of the N-order stationary 
points do survive the passage to the next order (N + 1), thus arriving for N — >■ 
oo at the exact RST spectrum in the spherically-symmetric approximation. 
The RST spectrum of order N = 3 is sketched in fig. 2 (at the end of the 
main text). 
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II Relativistic Two-Particle Theory 

Although the subsequent discussion refers exclusively to the non-relativistic 
limit of RST, it is important to realize that this theory is of truly relativistic 
character. Therefore the non-relativistic results displayed below are to be 
understood in the sense that they represent merely the non-relativistic limit 
forms of the corresponding relativistic outcomes which will however not be 
presented in the present paper. Nevertheless it is very instructive, for a bet- 
ter understanding of these non-relativistic results, to start from the original 
relativistic form of the RST dynamics and to glimpse at the emergence of its 
non-relativistic approximation which will then afterwards be inspected and 
treated in great detail. 

1. Hamilton-Lagrange Action Principle 

Similarly to most of the modern particle theories, the RST dynamics can 
also be based upon an action principle [T3] 

^W^RST = , (II.l) 

where the action integral PFrst refers to some Lagrangean /^rstI^,-^/^] 

WnsT = jd''xCnsT['^,A^] . (11.2) 

This Lagrangean appears as the sum of a matter part £d and a gauge field 
part £g 

Crst[^,A,] = Cj,m + Cg[A^]. (11.3) 
The matter part £d is given in terms of the wave function \& by [3] 

iJl/C — — — 

Cj^[\ij] = — [^r^' (p^,^) - (D^,^) r^^] - ^TWc^^ , (11.4) 
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and similarly the gauge field part Cq reads in terms of the bundle curva- 
ture {''field strength'' ) 

he 

Cg[A,] = K^pF\yF^^'' (II.5) 

he 

The curvature components do refer here to the decomposition of the 

bundle curvature jF^j, 

= - V,A,, + [A^, A] (11.6) 

with respect to a suitably chosen basis {tq} of the structure algebra 

This basis {tq,} builds up also the metric Ka/3 for the Lie algebra bundle in 
the following way: 

Kaf3 = ci tr(r„) ■ ti^Tfs) + Ca tr(ra ■ r^) . (II. 8) 

Here the coefficients ci and C2 are ordinary constants over space-time so 
that Kai3 is actually covariantly constant 

D^K^(3 = . (II.9) 

The gauge forces among the particles are incorporated into this formalism 
via the principle of minimal coupling, i.e. the covariant derivative of the wave 
function \1/ in the matter Lagrangean £d ( III.4P is defined by 

V^^ = d^^+A^^ . (11.10) 
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Here the bundle connection {^^ gauge potential") can also be decomposed 
with respect to the Lie algebra basis {t^}, quite analogously to its curva- 
ture j-^, (unD 

A = ' (11-11) 

i.e. for a two particle system whose gauge algebra turns out to be four- 
dimensional (a = 1 ... 4) with T3 = X, r4 = x- 

A = + ^%^2 + B^x - B^x (11. 12a) 

= F\,n + F\,T2 + G,,,x - *G^.uX ■ (n.l2b) 

These decompositions of the gauge potential and field strength demon- 
strate that the totality of gauge forces in RST subdivide into two classes: 
(i) electromagnetic forces being described by the real- valued four-potentials 
A°- ^ and field strengths F°-^y (a = 1, 2) and (ii) exchange forces being de- 
scribed by the complex- valued four-potential B^. However the point with the 
exchange type of interaction is that this kind of force can be active exclu- 
sively among identical particles [3[ HI] while in the present paper the interest 
aims at positronium as a bound state of two different particles, i.e. positron 
and electron. Therefore, for such a two-particle system, the exchange poten- 
tial and its field strength G^^ must be put to zero (5^ = 0, G^^ = 0) so 
that the gauge fields (lII.12aP - (]II.12bp are reduced to their Abelian projec- 
tions (a = 1,2) 

A\Ta (11.13a) 
J'f.u F\,Ta , (n.l3b) 
with the usual link of potentials and field strengths 

F\, = V,A\-V,A\. (11.14) 

14 



A similar simplification must occur also for the matter fields because the 
material particles (i.e. electron and positron) cannot feel the exchange forces. 
More precisely, since the two-particle wave function is the Whitney sum 
of the two single-particle wave functions ipai^ = 1,2), i.e. 

^(x) = ^i{x)®^2{x) , (11.15) 

the covariant derivative of \E' appears as 

V^^ = D^^Ij,®D^iIj2 (11.16) 

with 

D^^i = d^^i-tA\^, (II. 17a) 

D^^2 = ^,^2-^A\^2 . (n.l7b) 

Finally, it should be self-evident that the missing of the exchange forces 
does also simplify the fibre metric Kaj3 fill. 81) since this missing reduces 
the (4 X 4)-matrix Ka/B to the following (2 x 2)-matrix Kab 



{K^p} {Kab} = \ , (11.18) 




which obviously then simplifies the gauge field Lagrangean £g (III-5P to 



Cg[A,] ^ --^F\,F"^^ . (11.19) 
87rQ;s 
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2. Field Equations and Conservation Laws 

Once the kinematical setting for the action integral PVrst (HI .21) is estab- 
hshed, the dynamical equations of the theory do appear as the corresponding 
Euler-Lagrange equations; namely, the variation flll.ip with respect to the 
matter field yields the Dirac equation (DE) 

ihT'^V^^ = Mc"^ , (11.20) 

or, resp., in component form for the presently considered two-particle sys- 
tems [7] 

ihYD^ipi = -Mp cipi (II.21a) 

ih-f^D^^2 = MeCilj2 (II.21b) 

where Mp (Me ) is the rest mass of the positively (negatively) charged particle. 
(For the positronium system, both masses are identical: Mp = Me = M). 
Quite analogously, the variation flll.ip with respect to the bundle connec- 
tion yields the (non-Abelian) Maxwell equations 

P^^^, = -47rzas J. , (11.22) 

or, resp., in component form 

= 47rasJ^ (11-23) 
(a = 1...4) . 

But the point with these general gauge field equations is now that they 
simplify to the Abelian form 

V^F%, = 4na,j\ (11.24) 
(a = 1,2) , 
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because the absence of the exchange interactions for non-identical parti- 
cles (-W = 0) eliminates the non-Abelian part of the gauge field equations, 
see ref.s [3l [7]. 

The coupled set of matter and gauge field equations (IIL20p and (111.220 
is not yet a closed system because it is necessary to specify the Maxwell 
currents j"^ as the components of the current operator J'^ (]II.22|) 



^ ^J'^M^" = ^ (-^'^^1 + 3\j2 + c/^X - *g^x) , (n.25) 

namely just in terms of the wave function \E'. This requirement, however, can 
be satisfied by first linking the Maxwell currents j""^ to the RST currents ja^ 
by means of the fibre metric Kaf3 fill. 81) 

= K'^^jp, (II.26a) 
3a, = K^pj\ (II.26b) 

The point with this construction is namely that the RST currents Jq,^ are 
defined in terms of the wave function \I' through 

j,^ = ^T^r,^ , (11.27) 

and therefore these currents ja^i do automatically obey the source equations 

D^je,^ = (11.28) 

which therefore appear as an immediate implication of the matter dynamics. 
On the other hand, the covariant constancy of the fibre metric (111.90 
guarantees that also the Maxwell equations j"^ do obey a source equation, 
i.e. 

D^j\ ^ , (11.29) 
17 



which is nothing else than the component form of the abstract version 



V^'J^ = . (11.30) 

This identity, however, is itself an immediate implication of the gauge field 
equations flll.22|) so that the matter and gauge field dynamics are perfectly 



compatible with each other. 

Moreover, the absence of the exchange interactions for the positronium 
system reduces also the source equation (]II.29|) to its Abelian truncation 
(a = 1,2) 

V^3\ = (11.31) 

which itself is again an indispensable consistency condition for the Abelian 
simplification flll.24p of the general Maxwell equation (1II.22|) . Resorting here 



to the usual Lorentz gauge condition for the electromagnetic potentials A""^ 
(I11.13al) 

V^A\ = , (11.32) 

with observation of the curl relation (111.14p . converts those Abelian Maxwell 
equations (111.240 to the d'Alembert equations for the potentials A'^^ 

d^d^A\ = 47ca,j\ . (11.33) 

Thus the RST dynamics is revealed as a consistently coupled set of matter 
and gauge field equations, where the coupling of the gauge fields (a = 1, 2) 
to the single-particle wave functions ipa can be made somewhat more obvious 
by introducing the Dirac currents ka^ in the usual way 

ka/M = i'aliM^a ■ (11.34) 
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All three types of currents {j"^, ja^t, ^a/^} are then related to each other by 

= -J2M = (II.35a) 

j\ = -Ji, = -h, (n.35b) 

and appear as being correctly matched to the positive charge of the first 
particle (a = 1, positron) and to the negative charge carried by the second 
particle (a = 2, electron). 

The conservation laws, such as those for the electric charges (111.3ip . are 



relevant also in another respect: For a closed system with total energy- 
momentum tensor '•"'"■'T^i/ one expects the validity of the conservation law 

= . (11.36) 



For a field theory which is based upon an action principle flll.ip . such an 



energy-momentum conservation law like 011.36P must emerge via the well- 
known Noether theorems [13]; or alternatively it must be possible to deduce 
this also directly from the underlying dynamics (]11.20p and (]I1.22p [15j . Since 
the RST field system is built up by a matter part and a gauge field part, it 
should appear natural that the total tensor ^"'"''T^,y turns out as the sum of 
the matter density (^°^T^i.) and the gauge field density ("-"^^T^i.), i.e. 

= (°)T,. + (G)T,. , (11.37) 

with the matter part being given by 

(°)t^, = ^ [^r^ (p,^) - (D,^) r^^ + fr, (d^^) - (p^^) r,^] 

(11.38) 
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and the gauge field part by 



(^^T,. = ^ K^p {f\,F^' - \g,.F\^F^'^'^ . (11.39) 

By means of tlie dynamical equations for tlie matter fields 011.201) and gauge 
fields flll.22p one can now show that the conservation law flll.36p is actually 



obeyed because the sources of the partial densities *^°^T^,y and '^'^''T^y do just 
cancel 

(D)r^^ = {o)rj.^^ ^ flcF%,j^'' . (11.40) 

As we shall readily see, the existence of such a total density ^'^^T^y is impor- 
tant for the subsequent construction of the principle of minimal energy which 
will be based upon the total energy density ^'^^ToQ^r) as the time-component 
of (^)T,.. 



3. Stationary Bound States 

When the positronium system has settled down to its groundstate, the 
corresponding field configurations must be expected to be stationary or even 
time-independent. For the wave functions ipair^t) this implies the following 
product form (a = 1, 2) 

/ M \ 

V'a(r, t) = exp -t t ■ ^,(f) , (11.41) 



while the currents and gauge potentials will be found to be time-independent: 

ka,{r, t) {(")A;o (f) , -fc,(f)} (II.42a) 
A%if,t) {("Uo(f); -l.(f)} . (II.42b) 
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Accordingly, the d'Alembert equations (1II.33P are reduced to the ordinary 
Poisson equations 

A(^ko(r) = -47ras (^^fco(f) (IL43a) 

A(2ko(f) = 47ras ^^^A;o(r) (II.43b) 

AAi (f) = -Ana, ki (f) (II.43c) 

Ala if) = Aira, (r) (II.43d) 

whose formal solutions are given by 



li/2(f) = ±«s y rfV-'^;^ . (II.44b) 

Similarly, the time-independent mass eigenvalue equations for the spatial 
parts of the wave functions ipai'f^) (111.4ip may be deduced directly from the 
matter equations (111.21al) - (111.21bl) . or they can be deduced alternatively by 
substituting the stationary ansatz (1II.41I) into the action principle ( 111. II) - (111. 2^ 
and carrying through the variational procedure. The result is then found to 
be of the following form 

^7• V7/'i(r) + (%(r)7Vi(^) - Mr) ' iMr) = - (^7o + ^ j Mr) 



i^'V^r) + ^^%{r)l^i'2{r) - Ai{r)' i^r) = - ( ^^7o - ) Mr) ■ 



(I1.45a) 
(11.45b) 

For obtaining here the solutions ipair), it is very convenient to conceive the 
Dirac four-spinors ipa{r) as the (local) direct sum of the Pauli two-spinors 
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liJaif) = ('^V+lr") © ^-^V-lrl , (11.46) 

and for the subsequent treatment of the positronium groundstate one decom- 
poses these Pauh spinors with respect to a certain basis system {uj^\uj^{^'^} 
as 

W^+(f) = (r sint9)-5 (^('^)^+(r, d) ■ J^^ + ^''^S+{r, d) ■ 4"^) (II.47a) 
W(^„(f) = -i(rsin?9)-^ (^(")i?„(r, t9) ■ + (")^_(r,^9) ■ cu[^) . (II.47b) 

(For the details see ref. [7j). Clearly by such an arrangement, the original 
mass eigenvalue equations flIL45al) - flII.45bl) transcribe to the corresponding 
equations for the (real- valued) wave amplitudes ^^^R±, ^°'^S±] e.g. for the first 
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particle (a = 1, positron): 



dr r 

Mp - Ml 



h 

(IL48a) 



dr r d'd 

_ Mp - Ml 



c ■ (i)^_ 



h 

(II.48b) 

1 5(r(i)i?_) 1 



(^Uo ■ (1)^+ + (2)^^ (sin^9 ■ (i)^_ - cos^9 ■ 



r dr r d'd 

Mp +Mi 



(II.48c) 

r ov V 



r 9r 



Mp +Mi ^ (i)g 



(II.48d) 



Here it can be shown that for real-valued wave amplitudes ('')^±, ('^)^± 
the Dirac three-currents ka{r) 0II.42aP encircle the z-axis {d = 0, vr) of the 
spherical polar coordinates (r, -(9,0), i.e. 

ka{^ = ^^^k^{r,^)-e^ , (11.49) 

and therefore it is self-suggesting to adopt a similar symmetry for the mag- 
netic potentials A„(r) flll.42bl) : 

l,(f) = (%(r,^)-e;. (11.50) 
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The magnetic Poisson equations (111.43cp - (111.43dp are then transcribed to the 
azimuthal components ^""^A^ of these vector potentials Aa{r) as 

(1/2)4 

A(^/% - = t4™. . (11.51) 

sm V 

4. Non-Relativistic Approximation 

For the subsequent numerical demonstrations, we will restrict ourselves 
to the non-relativistic approximation; and for this purpose one therefore has 
to elaborate now the non-relativistic limit form of both the electrostatic Pois- 
son equations (111.43ap - (111.43bp and the mass eigenvalue equations (lII.48aP - 
flIL48dp . For the case of the Poisson equations, this problem is much simpler 
than for the eigenvalue equations, because the general form of the Pois- 
son equations survives the non-relativistic limit so that it becomes merely 
necessary to specify the non-relativistic approximation of the Dirac densi- 
ties ^""^ko^r) (111.42ap . But here one can easily show, see ref. [2], that the 
desired non-relativistic form of these charge densities is obtained by simply 
neglecting the "negative" Pauli components '■'^V-(^) (1II-47bp and by addi- 



tionally omitting alternatively one of the two spin directions. This means 
that the spin of any particle points to a definite z-direction, either the pos- 
itive one (^"-'>S+ = 0) or the negative one = 0); and thus both spin 
orientations become decoupled in the sense that the wave amplitudes do 
obey separate non-relativistic eigenvalue equations, e.g. for the first particle 
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(with positive z-component of spin, say): 



2Mn 



ld_ 

r dr 



dr 



(11.52) 



£"3(1) ■ 



(Es(i) = -(Mp +Mi)c^) 



or analogously for the second particle (with negative z-component of spin, 
say) 



^2 



2M, 



r dr 



dr 



J.2 Q^2 



(11.53) 



-^S(2) ■ ^^^'S^-i 



(ES(2) = (M2-Me)c2) . 

However, observe here that these non-relativsitic energy eigenvalue equations 
are not of the usual Schrodinger form! This is due to the fact that the 
eigenvalue of the total angular momentum (J^) is zero (see the discussion 
of this effect in ref. 0). But on the other hand the non-relativistic energy 
eigenvalues -Es(a) are linked here as usual to the difference of rest masses Mp/e 
and mass eigenvalues Ma as shown above. 

Concerning now the non-relativistic approximation of the electrostatic 
Poisson equations (lII.43aP - (lII.43bp one finds for the Dirac densities *^"^A;o {r) [2] 



^ >ko{r,§) ~ — — - 

47rr sm v 

(2), t (')^+2(r,^) 
Anr sm v 



(II.54a) 



(II.54b) 
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provided one adopts again the first (second) spin pointing to the positive 
(negative) z-direction. Thus the electric Poisson equations (III.43aP - (III.43bp 

appear as 

A (%(r, ^) = -a. . ^'^^^'("-f^ (n.55a) 

A(%(r,^) = as-^^^^^±^. (II.55b) 

r sm v 

In contrast to their electrostatic counterparts, the magneto static Poisson 
equations (lII.43cP -( ]II.43dp are disregarded because the magnetic effects ap- 
pear to be of the same order of magnitude as the relativistic effects, which 
of course are to be neglected for the present non-relativistic approximation. 
This is also the reason why the magnetic potentials ^°^A^{r^-d) have been 
omitted for the non-relativistic approximation of the mass eigenvalue equa- 
tions, whereas the original relativistic form (lII.48aP - (lII.48dp must of course 
contain the magnetic potentials ^°''^A^{r^'d) which, e.g., describe the spin-spin 
interactions for the groundstate. 

However, the crucial point with the residual electric Poisson equations 
flll.55al) - flll.55bp is now that they obviously will yield angular dependent po- 
tentials (")y4o(r, t?), even if one adopts the wave amphtudes '•^^-R+, '■^^5+ to be 
spherically symmetric! Clearly, these angular dependent potentials ''"^Ao(r, d) 
as solutions of the Poisson equations (]II.55aP - (lII.55bp must then entail their 
angular dependence also on the wave amplitudes ^^^R+i ^"^^Sj^ as solutions of 
the (non-relativistic) energy eigenvalue equations (]I1.52p - (111.53p . But this 
angular dependence of the RST eigenvalue equations will be neglected sub- 
sequently by resorting to the ^^spherically symmetric approximation'^ [?]. 
Namely, by this approximation do our non-relativistic RST predictions stand 
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on the same footing as the conventional treatment of positronium which re- 
hes upon the spherical symmetry of the Coulomb potential for the internal 
two-body problem (see the corresponding remarks in the Introduction) . Here, 
the neglection of the anisotropic effect means for the (non-relativistic) energy 
eigenvalue equations flll.52p - flll.53p that one simply puts to zero the angu- 
lar derivatives of the wave amplitudes ^^^R+ and '•^^5'+; and furthermore one 
adopts also the spherical symmetry for the electric potentials which then 
recasts those eigenvalue equations to the following truncated form: 

(II.56a) 



2Mp r dr \ dr 



(II.56b) 

Observe here that, in contrast to the situation with these spherically symmet- 
ric eigenvalue equations, the corresponding symmetric approximations t"ly4o(r) 
of the electric potentials ^"•'74o(r, •§) cannot be obtained by simply omitting the 
angular derivatives in their Poisson equations flII.55aP -f jII.55bp . but the de- 
sired spherically symmetric Poisson equations must be deduced directly from 
the principle of minimal energy which is now to be considered in greater 
detail. 
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Ill Principle of Minimal Energy 

As pleasant as the existence of the action principle OII.ip for RST may 
appear, the action integral ( 111.21) itself is not equipped with an immediate 
physical meaning. In this context, it seems highly desirable to have some vari- 
ational principle which does refer to the total field energy being stored in 
the RST field configurations. If this goal could be attained, both the mass 
eigenvalue equations and the Poisson equations would not only appear as the 
(stationary) Euler-Lagrange equations due to the Hamilton-Lagrange action 
principle (]II.ip - (lII.2p but this eigenvalue problem would then appear also in 



form of the variational equations due to that total energy functional Erp. 
Possibly, the variations of the functional i?T must be restricted by some con- 
straints. Fortunately, it has already been demonstrated that the construction 
of such an energy functional is possible, for both the relativistic and the 
non-relativistic situation [31 [7] . 

1. Relativistic Construction 

Naturally, one will assume that the wanted energy functional (-Et) will be 
related in one way or the other to the total energy density E't as the spatial 
integral over the total energy density '•'^^Too(r): 

Et = j d^f (T)Too(f) = + (Ill.la) 
E^ = j rf¥(D)Too(f) (Ill.lb) 
Eg = j (iV(«)Too(r-0 . (III.lc) 

Clearly, since the total energy-momentum density '■^^T^jy fill. 371) appears as 
the sum of the matter part '•^•'T^jy flll.38p and of the gauge part ^'^'^T^i^ fill. 390 
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the total energy E-r is the sum of the matter energy E^) and the gauge field 
energy Eq. 

Furthermore, since the two-particle wave-function \l'(r) in the matter den- 
sity (°)T^^ fllL38D is the Whitney sum flll.lSp of the two one-particle func- 
tions ipai'T^, the matter energy Ed fllll.lbp is revealed as the sum of the 
individual energies £'D(a) of both particles: 

Eb = Ed (I) + Ed (2) , (III. 2) 

with the one-particle contributions -Ed (a) being given by [3] 

Ed(i) = Zfi) ■ Mp + 2Tki„(i) + Mt^c" (III.3a) 

Ed(2) = ^2) ■ Me + 2Tki„(2) + Mif ^c^ . (IIL3b) 

Surely, this is a very plausible result because it says that the material ener- 
gies -Ed (a) of both particles consist of: 

(i) the rest mass energies Mp/eC^, to be corrected by the mass renormal- 
ization factors Z(a) 

Zl) = j d?ri,,{r)i;,{r) (III.4) 

(ii) the proper kinetic energies Tkin(a) 

Tkin(i) = '-he j rf¥V'i(r)7- (III.5a) 
Tkin(2) = -\nc j t/¥7/>2(r)7- W2(r1 (111.5b) 



(for the factor of two in front of the kinetic energies Tkin(a) ( IHI.SaP 
(lIlLSbl) see ref. [16]) 
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(iii) the mass equivalents Mj™jC^ of the magnetic interaction energy 

j^^M^2 _ J d3^X2{r) ' ki{r) (III.6a) 

M^^'^c^ = hc j d^rAi{r) • k2{r) . (III.6b) 

Observe here that, in contrast to the magnetic case, the mass equiva- 
lents MjJjjC^ of the electric interaction energy 

M^'^c^ = hc j rf¥(%(f) • (i)fco(f) (III.7a) 

Mj^^c^ = -he [ rf¥(%(f) ■ (%(f) (III.7b) 



are not part of the matter energies -Ed (a) cind thus there is here no perfect 
symmetry between the electric and magnetic interactions. This symmetry, 
however, is manifestly present in the gauge field Eg fllll.lcp which does ap- 



pear here as the interaction energy of the real gauge field modes, i.e. 

with the electric and magnetic parts being given by 

^(e) ^ _nc_ I ^3_.^(i^ ^ V(%(f) (111.9a) 
47ras J 

.(m) _ f 3 



E^^'"' = / d'fHi'H2 (111.9b) 



(#„(f) = V X l,(f) 



But though, with these arrangements, the total field energy Et (1111. laD 
is well in agreement with the conventional field-theoretic conceptions, this 
functional E^ cannot serve directly for the wanted variational principle. The 
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reason is that the corresponding variational equations (due to SE^ = 0) do 
not coincide with the coupled set of Poisson and mass eigenvalue equations 
flII.43aD - flII.43dp and flII.45aD - flII.45bp ! The reason for this is mainly that the 



total energy functional i?T flllLlaD - flllLlcP contains no coupling of the mat- 
ter fields ^a{r) to the electric gauge potentials ^"'^AQ{r), while such a coupling 
is present in both the mass eigenvalue equations fill. 45aP - f ill. 45bp and the 



Poisson equations flII.43aP - flII.43dp . Consequently, it becomes now necessary 
to introduce the missing couplings in such a way that, on the one hand, 
the numerical value of the original functional i?T upon the solutions of the 
RST eigenvalue problem is left unchanged but, on the other hand, the cou- 
pled eigenvalue and Poisson equations do actually emerge as the variational 
equations due to the modified functional E^. This problem of modifying the 
original E^ to i?T under preservation of its numerical value can obviously 
be solved by imposing upon the variational process certain constraints which 
are automatically obeyed by all the solutions of the RST eigenvalue problem. 
Here, the first constraint refers to the fact that the mass eigenvalue equations 
f lll.45ap - flll.45bp for fixed potentials ^"'^AQ{f) , Aa{f) are formally linear which 



admits to normalize the solutions ipai'f^ to unity (a = 1,2): 

(i¥(%(f) = [ d¥VI(r)^„(f) = 1 . (III.IO) 



Thus the first constraint is that of wave function normalization and reads as 
follows: 

iVD{a) = J d^fM^^'H^ -1 = . (III.ll) 

The second constraint to be imposed refers to the Poisson equations 
f lll.43ap - flll.43dp . Multiplying through these equations by the corresponding 
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potentials and integrating by parts over whole 3-space (under use of Gaufi' 
integral theorem) yields the following Poisson identities [3l [7] 



(e) 



G{2) 



N, 



(m 



G(2) 



d^r 



d^r 



V (^ko(f) • V(%(rO + Anas ^%{^ ■ ^%(r) 



= 



= 



V X y4'i(f) ) • (V X A2{r) ) + Anas Ai{r) • ^(^ 



(III. 12a) 

[III.12b) 


(III. 12c) 



V X Aiir)) • (V X Ai{r)) - Ana^ Ai^r) • ki{r) 



. 
(III.12d) 

These two kinds of constraints (i.e. wave function normalization and Pois- 
son identities) must now be implemented into the naive variational princi- 
ple {6Et = 0) by means of the method of Lagrangean multipliers so that 
one arrives at the modified variational principle ( principle of minimal 
energy) : 

SEt = (III. 13) 

with the modified energy functional i?T being defined as the sum of the 
original functional f llll.lap and of the constraints flIII.lip - flIII.12dp : 



(e) _ ^(m) 
R 



Et = ■ Mp + Zf2) ■ Me C^) + 2 (Tkin(l) + Tkin(2)) + ( " ^ 



a=l 



;iii.i4) 
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Here the magnetic mass equivalents Mj™jC^ occurring in the original func- 
tional fllll.lap . namely via the matter energies -Ed (a) (!III-3ap - fllll.3bp . 
have been eliminated in favour of the magnetic field energy E^^^ fUTTObl) by 
reference to the magnetic Poisson identities (lIIL12cP - (lIII.12dl) which read in 
physical terms 

= -Mi'^^c^ = -M'lfc^ . (III. 15) 

But here it is now a nice exercise to convince oneself of the pleasant 
fact that both the relativistic mass eigenvalue equations (lII.45aP - (]II.45bp and 
the Poisson equations flll.43al) - flll.43dp actually turn out as the variational 
equations of the functional (1111.14^ . provided one fixes the Lagrangean mul- 
tipliers Ad (a) and A^^^'' as follows: 

Ad(i) = Mic^ (III. 16a) 

Ad(2) = -Msc^ (111.16b) 



)(2) 

iH , _ \(m) 



he 



(a= 1,2) . 



Concerning the general form of our final result E't f lIII.14p it is interesting 



to observe that this functional is built up by two different kinds of contribu- 
tions, namely the physical terms (first line) and the constraint terms (second 
line). Thus, when one wishes to resort to some trial configurations for the pur- 
pose of minimalizing the functional Et, it is advantageous to work with those 
ansatz functions which obey both types of constraints (lIII.lip and (1III.12aP - 
(IIII.12dl) . In this case, the second line on the right-hand side of equation 
(IIII.14P is zero and the functional E^ becomes built up exclusively by the 
physical terms (first line). Both cases (i.e. with vanishing and non- vanishing 
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constraint terms) have been tested in the preceding papers 0[ Ej and it has 
been found that the theoretical predictions are in better coincidence with the 
experimental data when the constraint terms are made zero by an appropiate 
choice of the trial functions (see ref. |3j). But if one should prefer (from some 
reason) a variational ansatz which implies non-vanishing constraint terms, 
then the non-trivial contribution of the latter terms ensures that the energy 
prediction is kept in the vicinity of the proper energy value see ref. [7]. 
Subsequently we will use trial configurations with vanishing constraint terms 
and will then study the associated (albeit non-relativistic) energy predictions 
for positronium. 

2. Non-Relativistic Approximation 

The non-relativistic two-body problem is perfectly understood by the con- 
ventional theory, both classically and quantum-mechanically [8]; and there- 
fore any new relativistic theory must not be admitted to come into conflict 
with the known (non-relativistic) results of the conventional theory. For the 
present situation this means that the non-relativistic positronium spectrum, 
as predicted by RST, should agree with the corresponding conventional pre- 
dictions (at least up to the relativistic corrections). Thus it becomes neces- 
sary to look for the non-relativistic approximation (e!^^) of the relativistic 
energy functional E^ f lIII.14p . so that the non-relativistic positronium spec- 
trum can be estimated by means of appropriate (but non-relativistic) trial 
functions. 

The non-relativistic form e!^^ of the energy functional £'t is, however, 
easily deduced by simply looking for the non-relativistic forms of the indi- 
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vidual contributions. Thus, omitting as usual the rest masses for the non- 
relativistic approach, the first two bracket terms on the right-hand side of 
equation flIII.14|) yield the ordinary kinetic energies -Ekin(a) of the two parti- 
cles 

{■^(1) ' ^ + ^(2) ■ C^) + 2 (rkm(l) + ^km(2)) ^ -£'kin(l) + -£'kin(2) ; 

(III.17) 

with the non-relativistic energies being given by [2] 



km(l) 



kin (2) 



4M„ 





I dr 


2 

) +^ 




(111.18a) 




\ or 1 


2 




(111.18b) 



{cff = rdrM) . 

Next, neglecting the magnetic interactions for the present non-relativistic 
limit, the general form of the electrostatic field energy E^ (lIII.Qal) remains 
the same for the non-relativistic limit; but the electric counterpart of the 
magnetic Poisson identities, i.e. 



E, 



(III.19) 



receives some change concerning the electric mass equivalents My\^c^ flIII.7aP 



flIII.7bp . The point here is that the original relativistic form of the Dirac 
densities ^^^/cq (r) 

(a)6 2 ^ (a)^_2 ^ {a)g 2 ^ {a)gj 



(rl 



(III.20) 



47rr sin d 

becomes truncated to the first contribution (~ ^"•'i?+^) as shown by equations 
( ]II.54aP -( lII.54bl) . Indeed multiplying through the non-relativistic Poisson 
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equations (111.55ap - (111.55bp and integrating by use of GauB' theorem yields 
the non-relativistic form of the electric Poisson identities piI.12aP - (IIII.12bp 

as 



^G(l) =^ ^G(l) 



rf¥ V(%(f)- V(%(r) + 



a. 



= 



^^G{2) 



^G(2) = I d'r fv(%(f).V(%(f)- 



r sin-i? 

(IIL21a) 



r sin?? 



. 



[III.21b) 



Thus the final form of the desired non-relativistic energy functional e'^'^ 
is the following 

2 

E!^^ = -Ekm(l) + -^^11(2) + -^R^ + ^ \^S{a) " !«) + -^G^a) ' ^G 



a=l 



(III.22) 



Here the non-relativistic kinetic energies -Ekin(a) have been specified by equa- 
tions flIII.18aP - flIII.18b[) : the non-relativistic electric field energy E^^ pre- 
serves its relativistic form fllll.Qal) but is now related to the electric mass 
equivalents MjJjjC^ flIII.7ap - flIII.7bp in their non-relativistic form, i.e. 



Anr sin 



Mi^^c" ^ -he J d¥(%(f) • 



Attt sin -(9 



(III.23a) 
(III.23b) 



and the electric Poisson constraints ^Vq^q) have to be used in their non- 
relativistic form N§\a) flIII.21aD - flIII.21bD . Finally, the non-relativistic form 
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of the normalization conditions (lIII.lip is now 



A^D(i) ^ j d'r -1=0 (IIL24a) 

r (2)0 2(^\ 

iVD(2) < 2) = / d'f -1 = 0, (III.24b) 

w J 47rr sm v 

with the relativistic Lagrangean multiphers Ad (a) adopting their non-relativistic 
form As (a) as the Schrodinger energy eigenvalues -Es(a) ( III.52p - (]lI.53P 



2/ 



Ad(i) = Mic^ ^ As(i) = -^s(i) = Mp + Mic^ (111.25a) 
Ad(2) = -M^c^ As(2) = -^s{2) = Me - MaC^ . (111.25b) 

But with all these arrangements it becomes now a straightforward mat- 
ter to convince oneself that the non-relativistic eigenvalue equations ( 111.521) - 
( 111.531) together with the non-relativistic Poisson equations (111.55a|) -( ]11.55b[) 
actually turn out as the variational equations {6e!^^ = 0) due to the non- 
relativistic energy functional ^r^*''*! However despite of this pleasant result, 
we will not use this functional e!^^ (1111. 22p in its present form for the subse- 
quent calculation of the non-relativistic positronium spectrum but will first 
resort to a further simplification. 

3. Spherically Symmetric Approximation 



The wave amplitudes ^"'^R±, ^"'^S± (]11.47aP - (]11.47bp are assumed to depend 



not only upon the radial variable (r) but also upon the angle i) of the spher- 
ical polar coordinates. This angular dependence must be taken into account 
also for the positronium groundstate (as well as for the excited states of the 
same symmetry) despite the fact that the groundstate is mostly the state 
of highest symmetry, i.e. the spherical symmetry for the present situation. 
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It is true, the non-relativistic eigenvalue equations (111.52p - (111.53p would ad- 
mit spherically symmetric solutions ^^''-R+(r), *^^^5-|.(r) provided the electric 
potentials ^"•^Ao{r,'d) are also symmetric: '•"^y4o(r, -(9) =^ t"Uo(r); but even for 
SO (3) symmetric amplitudes ^^^R+{r) and '^^o+(r) the potentials ("Mo(0 as 
solutions of the non-relativistic Poisson equations flll.55ap - flll.55bl) can then 
obviously not be spherically symmetric. But because of these anisotropic 
potentials, the solutions and of the eigenvalue equations will then 
also be found to miss the spherical symmetry and this contradicts the origi- 
nal assumption for the wave amplitudes. This anisotropic feature of the RST 
groundstate solutions is in contrast to the SO (3) symmetry of the residual 
internal Coulomb force problem in the conventional theory after the centre- 
of-mass motion has been separated off [6j. Thus, in order to compare the 
RST predictions for the non-relativistic positronium states to the correspond- 
ing conventional predictions on the same level of symmetry, one will restrict 
the RST trial functions to just this spherical symmetry. Furthermore, one 
assumes that both particles do occupy physically equivalent states, i.e. the 
wave functions and potentials due to both particles must be identical. This 
yields the following identifications 

(^)^+(r, = (^)5+(r, ^) = R{r) (III.26a) 

Wfco(r, ^) ^ (%(r, ^) = -^Mi = (iii.26b) 

Arcr sm v Anr sm ly 

(^Uo(r, ^) = -^^^Ao{r, ^) = [PUo(r) . (III.26c) 
Accordingly, the two non-relativistic eigenvalue equations pi.52p and flll.53p 
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become contracted to a single one 



2M 



Id dR{r) 
r dr \ dr 



- he [PUo(r) ■ R{r) = Es ■ R{r) . (III.27) 



On the other hand, the deduction of the spherically-symmetric form of 
the two associated Poisson equations flll.55a|) - flll.55bl) is not so evident and 
requires an extra argument: For obtaining the S0(3) symmetric form of the 
non-relativistic energy functional E^^^ (IIII.22D one observes first that by the 
present identifications (]III.26aP - (lIIL26bp the kinetic energies -Ekin(a) (Illl.lSaP - 
(jlll.lSbp become identical 

2 



km(l) 



E, 



kin (2) 



AM 



" , / dR(r) 
drr 



dr 



(III.28) 



Next, the spherically symmetric form of the electric field energy aj 
becomes by reference to the isotropic potential [PlAo(r) (]III.26cP 

2 



E, 



(e) 
R 



he 



a. 



drr 



2 / d^oir) 



s Jo 



dr 



(III.29) 



and finally the non-relativistic normalization conditions (]III.24aP - (]III.24bp 
contract to a single one for the residual wave amplitude R{r) flIII.26al) 



(1) - iV^ (2) - iVj3 - 



dr rR{r) 



1 = 



^111.30) 



Of course, both non-relativistic Poisson identities flIII.21aP - flIII.21bl) do also 
collapse to a single one: 



7Vq (1) - (2) - 1\q - 



. 



An I drr"^ 
'o 



d^ 
dr 



TT [PUo(r) ■ R(rf 
2 r 



(III.31) 
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Thus the spherically symmetric form of the non-relativistic functional e!^^ 



(IIII.22P is found to look as follows 



2 M Jo y dr J a, Jo \ dr j (111.32) 

+ 2As-<)+2Ag^-iv(°\ 



where the Lagrangean multipliers for the wave function normalization flIII.24aP 



( ]III.24b[) become identical, too, and equal then the common non-relativistic 
energy eigenvalue Ea, 

-^S(i) = -^S(2) ^ -^s = — -E's (III. 33) 

(the multipliers flIII.16cl) for the electric Poisson constraints are identical in 
any case (a^V)) = ^§^2) = Ag^)). 

Now it is again an easy exercise to convince oneself of the fact that the 
SO (3) invariant functional eI^^ ( ]III.32p actually is stationary upon the so- 



lutions R{r) of the SO (3) invariant eigenvalue equation (1IIL27P as required. 
But the point here is that the desired Poisson equation for the rotationally 
invariant potentials tPly4o(r) is obtained now as the corresponding variational 
equation of the present energy functional flIII.321) : 

r'^ dr \ dr J 2 r 

This completes the non-relativistic, spherically symmetric RST picture of the 
positronium system. It consists of the eigenvalue equation for the common 
wave amplitude R{r) flIII.271) and the Poisson equation flIII.341) . with the 
energy functional given by flIII.321) . The pleasant effect with this result is 
now that the energy functional ^r^^'' (1111.320 is stationary upon the solutions 
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of the coupled matter-gauge field system (1III.27P and (1III.34P , so that one can 
now guess appropriate trial functions for the wave amplitude R{r) with ansatz 
parameters bk in order to extremalize the energy functional E^?^ flIII.32D just 
with respect to these trial parameters bk- 

4. Poisson Identities 

As mentioned above, the accuracy of the RST predictions for the total 
energy i?T become improved when one does not resort to two rather unre- 
lated ansatze for the wave amplitude R{r) and for the electric gauge poten- 
tial fp]ylo(r) but when one observes the link between both ansatze which is 
established by the electric Poisson identity Nq Hi11.31\} . or by (lIII.lQp . resp. 
Actually such an improvement can be attained by choosing some trial ansatz 
for the wave amplitude R{r) (as realistic as possible) and associating to it 
the corresponding solution ^^^Ao{r) of the Poisson equation flIII.341) . Subse- 
quently, we will elaborate this point in great detail by means of a special trial 
function R{r). 

Here, it suggests itself to try an exponential (~ exp(— r/r*)) times some 
other function {P{r), say), i.e. our general trial ansatz will look as follows 



Rir) = ^ —P{y) exp{-(3y) , (III.35) 

where r^, is some constant length parameter which mainly serves to obey the 
(non-relativistic) normalization condition (lIII.30p . The radial variable r of 
the spherical polar coordinates is rescaled to the dimensionless variable y 

y = — (111.36) 

Ob 



f'^B = , , o • • • Bohr radius^ , 
V Me"' / 
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and the decay constant 13 is to be considered as one of the ansatz param- 
eters hk- The other parameters are thought to be contained in the ansatz 
function P{y) which is required to not tend to infinity faster than exp(/5y) 
for y —>■ oo, so that the whole trial function R{r) safely tends to zero at 
spatial infinity {y oo). But clearly before one looks now for the associated 
potential tPUo(r) according to the Poisson equation, one will first satisfy the 
normalization condition (]III.30P which adopts the following shape 



dyyQ{y)e-^^y={^-^)' (III.37) 



where the constant is related to the original normalization parameter 
flIII.35p through 



y. = — (III.38) 

Ob 

and the function is the square of -P(?/), i.e. 

Q{y) = P{yf . (III.39) 

After the general form of the trial wave amplitude has thus been fixed, one 
next wishes to specify also the general form of the associated solutions [Ply4o(r) 
of the Poisson equation flIII.34p . For this purpose one first writes down the 
formal solution of this equation for the field strength [^'^(r) as 

WF(r) = ^^^^ = -^(l-- rdrrRirA . (III.40) 
dr \ 2 Jj. J 

Here one demands that the electric field strength (i.e. the gradient field of 
the potential ^P^ylo(r)) asymptotically adopts the usual Coulomb form due to 
one charge unit: 



lim WF(r) = . (III.41) 
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Clearly, this is easily seen to be a consequence of the normalization condi- 
tion fllll.301) . It is more convenient to recast the field strength ^^'^(r) to its 
dimensionless form 



[plF(r) = [i - (i + f(y)) e-H , (III.42) 

where the dimensionless function f{y) has evidently been introduced through 

poo 

f{y) = -1 + e^^y - dr' r'R{r'f . (111.43) 

^ Jr 

Observe here that, by virtue of just the normalization condition OIILSOp . 
the function f{y) vanishes at the origin {y = 0) so that the electric field 
strength ^^(r) (IIII.42D may adopt there a finite value (for an example of 
this see ref. [7]). One may therefore put 

fiy) = 2Py + Giy) (III.44) 

where G{y) must then tend to zero quadratically at the origin, i.e. 

Giy)=g2-y' + o{y^) (III.45) 
{92 = const.) . 



Namely, by this arrangement the value of the field strength WF(r) flIII.40p 
at the origin {y = 0) is then found as 



Wf(0) = -^(2/?2-^2) , (III.46) 
Ob 



whereas its general form is found to look as follows: 

WF(r) = • 1 (Piy) - G{y) ■ e'^^y) , (III.47) 

cb^ r V / 
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with the dimensionless function F{y) being given by 

F{y) = 1 - (1 + 2(3y) e'^^^ . (III.48) 

The field strength [plF(r) flIII.47p aheady displays its final form to be 
used subsequently for the computation of the field energy E^^^ fUTLM]) . and 



therefore it is advantageous to link the function G{y) directly to the ansatz 
function Q{y) (1III.39I1 . This, however, is easily achieved by substituting the 
field strength ^^'^(r) (1III.47I) back into the Poisson equation ( 1111.341) in order 
to deduce the desired differential equation for G{y) as 

-^^ + 2pGiy)=yQiy) (III.49) 

where the function Q{y) is a slight modification of the original Q{y) ( 1III.39P : 



2X2 



Q{y) = (^-J Q{y) - m . (iii.so) 

This new function Q{y) obeys the integral relation 

POO 

/ dyyQiy)e-^Py = 0, (III.51) 
Jo 

which of course is an immediate consequence of the former normalization con- 
dition ( 1III.37I) . For the subsequent demonstrations with appropriate ansatz 
functions Q{y), the associated constituent G{y) of the field strength [PlF(r) 
( ]III.47p is then easily calculated from the differential equation (1111.490 as 



POO 

G{y)=e^^y- / dy'y'Q{y')e-'^y' . (III.52) 
•>y 

However with respect to the Poisson identity, the most important feature 
of the field strength WF(r) ( 1111.47^ refers to its sum structure which implies 
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the splitting of the field energy E^^ flIII.291) into three contributions 



aB Jo \ y 



as Jo \ y J 

The first contribution (due to F{y)) essentially represents the field energy of 
the groundstate, whereas the second (quadratic) term and the third (quartic) 
term describe the infiuence of the excited states (see below). Observe also 
that all three integrals in flIII.53P assume finite values, namely since the field 



strength WF(r) 0III.47P remains finite at the origin {y = 0), cf. 0111.460 . 

With the field energy E^^ being analyzed now in sufficient detail, a similar 
investigation of the electric mass equivalents MjJjjC^ is necessary in order to 
elucidate the Poisson identities (1111.19^ more thoroughly. First observe here 
that, by virtue of the groundstate identifications (lIII.26aP - (lIII.26cp . both 
electric mass equivalents Myj^c^ (lIII.23ap - (]III.23bp become identical 

poo 

Mi'^c" = Mjf^c^ # M(^)c2 = --he / drr [PUo(r) ^% (r) . (III.54) 

2 Jo 

Here, the common reduced Dirac density ^''^ko (r) is given by the non-relativistic 
approximations together with the groundstate identifications as 

^%{r)=R{r)\ (III.55) 

and furthermore the electric interaction potential [p'y4o(r) is decomposed, 
analogously to its field strength MF(r) flIII.47p . into the groundstate con- 
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tribution A{y) and the excitation part A{y) as follows: 

[pUo(r) = ^ . 1 (A{y) - A{y) exp[-2/3y]) . (III.56) 



Both constituents A{y) and A{y) must of course be determined by identifying 
the derivative of the potential My4o(r) with its field strength WF(r) (1111.470 
which then yields the following two differential equations: 

A{y) - y^ = F{y) (III.57a) 

(1 + 2Py) A{y) - = G{y) . (III.57b) 

The first one (llll.57aP of these equations can be immediately solved be- 
cause the right-hand side F{y) is explicitly known, cf. fllll.48p . which yields 

A{y) = y r dy'^ = l- e-'^y . (III.58) 

Jy y'^ 

The second equation (]111.57bl) can of course not be solved explicitly since the 
right-hand side G{ii) will depend upon the selected trial ansatz, cf. fllll.52p : 
but for the practical purpose (see below) it is convenient to look for the 
formal solution A{y) 

A{y) = ye''' ■ r %e-'^y'G{y') . (III.59) 

Jy y 

Thus the excitation part A{y) of the potential fPly4o(r) (1111.56^ is non-singular 
at the origin {y = 0) and actually becomes zero: ^(0) = 0, according to the 
behavior of the function G{y) (1111. 45|) . This, however, then implies that the 
electric potential tP]Ao(r) (1111. 56|) itself is finite at the origin, i.e. 



[pUo(O) = 2(3-— (111.60) 

Ob 
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(for a sketch of this type of potential see ref. [7j). Amazingly enough this 
value at the origin does depend solely upon the ansatz parameter P, but 
not upon the other ansatz parameters bk which are contained in the ansatz 
function P{y) flIII.35p . Obviously, a pleasant property of these Coulomb- like 
potentials is that they yield a well-defined finite field energy m 
contrast to the situation with the exact Coulomb potential. Clearly, when the 
field energy E^^ would be infinite, such relations like the Poisson identities 
(IIII.19P would be meaningless! 

With the general shape of both the electric potentials A{y),A{y) and 
the trial function Q{y) being fixed now through the equations fllll.501) and 
f llll.58p - fllll.59p . the corresponding general form of the electric mass equiva- 
lent M^'^^c^ fllll.54p is found to look as follows: 

2 rco 

M^'^^c^ = -— . (2/3)' / dyA{y) ■ e'^^^/ 

O-B Jo 

dy Aiy)Qiy) ■ e"'/^^ + — ■ (2/3)' / dy i(y) ■ e"^^^ 

Ob Jo 

+ — dy A{y)Q{y) . e-'^y . 

«B Jo 

(III.61) 

Now according to the Poisson identities, this result for the mass equivalent 
must be exactly identical to the corresponding result for the field energy 
f lIII.53p ! But here it appears as a matter of course that the required identity 
of both results must already apply to the individual contributions themselves. 
Thus, considering first both groundstate contributions, one expects that the 
following identity holds: 

dy (^^^ = m' dy A{y)e-'Py . (III.62) 
47 
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flB Jo 



However, the validity of this relation is easily checked by simply substituting 
for the functions F{y) and A{y) those results obtained through equations 
flIII.48p and fllll.SSp which yields 

dy I ^1 = /3 (III.63) 



V y 

Next, the Poisson identity demands the identity of the quadratic terms 
in flITL53|) and fUTTeT]) . i.e. 



oo 



2. r'^e~'^y-F{y)G{y) = (2/3)^ H dy A{y).e~'^y- I dy A{y)Q{y).e 

Jo y Jo Jo 

(III.64) 

But here one substitutes for the function F{y) the derivative of A{y) from 
equation flIII.57aP and obtains through integrating by parts with regard of 



the differential equation for G{y) flIII.49P the following relation 



r ^ e-'^y ■ F{y)G{y) = - T dy e'^^^ ■ A{y)Q{y) . (III.65) 

Jo y Jo 



This recasts the required identity for the quadratic terms (1111.641) to the 
following form 

/ dy e-^^y ■ A{y)Q{y) = - (2(3)' / dy i(i/)e~^^^ (III.66) 
Jo Jo 

whose correctness is easily proven by use again of the link of A{y) to G{y) 
( ]III.57bp . Thus the quadratic terms of the field energy E^^ ( 1111531) are 
actually identical to those of its mass equivalent M^^h"^, fllll.GlD . 



Finally, it remains to prove also the identity of the quartic terms, i.e. 

dy =-J^"^dy A{y)Q{y) ■ e-'^y . (III.67) 
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The correctness of this equation is easily reahzed by substituting for the mod- 
ified function Q{y) from equation flIII.491) and integrating by parts with use of 
the differential equation for A{y) flIII.57bl) . This completes the demonstration 
of the electric Poisson identity in terms of the auxiliary functions Q{y), G{y) 
and A{y) which build up the interaction potential ^'^^AQ{r) flIII.56p and its field 
strength WF(r) (]III.47|) . Indeed for any trial function R{y), these three func- 
tions must be calculated explicitly in order to obtain the field energy E^^^ (or 
its mass equivalent M^^^c^) as part of the total energy e!^^ mrm . since 
the required trial functions R{r) may possibly look very complicated, the 
computation of the total energy e!^^ will also be somewhat intricate; and 
therefore the confidence in the results may be supported by testing explicitly 
the present identities flIII.65l flIII.66p ) and flIII.67p . This procedure is exem- 
plified now by a further specialization of the trial functions R{r) (]III.35p . 
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IV Hydrogen-like Wave Functions 

Surely, it will be very difficult (if not impossible) to find exact analytic 
solutions of the present (non-relativistic) eigenvalue problem ( 1111.271) plus 
(llll.34p . But nevertheless it may be possible to get the corresponding energy 
eigenvalue e!^^ very accurately if one succeeds to guess some trial function R 
as realistic as possible. Here it was observed that for the (non-relativistic) 
groundstate a simple exponential trial function R{r) is sufficient in order to 
reproduce exactly the groundstate energy Eq of the conventional Schrodinger 
theory, see below and ref. [3J. Incidentally such a simple exponential function 
represents the exact groundstate solution of the conventional theory which, 
by separating off the centre-of-mass motion, leads to the ordinary Coulomb 
force problem for the relative motion of both particles. Therefore it appears 
self-suggesting to select as the RST trial functions just that type of exact 
solutions for the conventional Coulomb force problem, i.e. the hydrogen-like 
wave functions. These are given by the product of some polynomial times 
an exponential function. In this sense, we specialize now the hitherto unde- 
termined ansatz function P{y) (1111. 35p down to a polynomial of degree N 

N 

p{y)^PM = J2^ny'' (iv.i) 

n=0 

where the coefficients 6„(n = 0, 1, 2 . . . N) work as the ansatz parameters (be- 
sides the parameter f3) and the variable y is the rescaled radial distance r 
(1111. 36p . Now, for inspecting more closely the groundstate and the excited 
states, it is very instructive for the practical use of the RST principle of min- 
imal energy to ffist convince oneself of the validity of the Poisson identities 
in terms of the polynomials Pn{y). 
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Turning here first to the modified square Q{y) (lIII.SOp of the polyno- 
mial Pn{y), one gets a polynomial Q2^{y) of degree 2N, i.e. 

N 

Q{y) Q2n{y) = Yl ^Jny^^"" - m' , (iv.2) 

m,n=0 

with the rescaled ansatz parameters b'^ being simply defined through 

&n = - ■ &n . (IV.3) 

y* 

Here, the constant y^ itself becomes a function of the ansatz parameters fe„ 
just by virtue of the normalization requirement flIII.37p 



m,n=0 ^ 

where the first coefficient 60 may be adopted to be unity (60 = 1), without 
loss of generality. Alternatively, the normalization condition pV.4p may be 
recast also to the following form 

.--(2,)^^- t yjj^^^^0^ (.V.5) 

n,m=0 ^ ' 

(n+m>l) 

which will be needed hereafter for the explicit calculation of the energy eigen- 
value eI^\ Observe also that, with the help of the latter form fllV.SP of the 



normalization condition, the modified square Q2n{y) f!IV.2p can be rewritten 
in the following form for N > 1 

N 

(y) (IV.6a) 

n,m=0 
(n+m>l) 

* - / - ^ (IV.6b) 
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which is best suited to verify the integral relation (lIII.Sip . namely through 
the somewhat trivial integral for any integer k: 

POO 

/ dyyqkiy)e-^^^ = 0. (IV.7) 

With a consistent variational ansatz being now at hand, one can proceed 
to compute the corresponding electrostatic objects in order to test the va- 
lidity of the Poisson identity just for that specific ansatz. First, turn to the 
function G{y) as a constituent of the electric field strength WF(r) (1111.470 and 
find from the constructive relation (11II.52P for the corresponding polynomial 
of degree 2N + 1: 

N 

Giy) G^n+M = -(1 + 2(3y) + b'Mm + n + 1)1 ■ ~g^„iy) (IV.8) 

m,n=0 

with the functions gmniy) being defined through 

9mn{y) = Y] 7 ■ — ■,Uoa\u+l ■ (^'^•^) 

^ {m + n + 1 — uy.{2py+^ 

But with the functions (52^+1 (HV-Sp and Q2n(|/) flIV.6ap - fllV.6bp being explic- 
ity known, one can check now the validity of the first integral identity f]111.65P 
and one actually finds its polynomial realization in the following form: 

r ^e-'Py ■ F{y)G,n+i{y) = - Hdye-'^y ■ Aiy)QMy) 

Jo y Jo 

= P 'Z^ PmPn{m + ny.Kmn ■ 

n,m=0 
(n+m>l) 

Here for the sake of simplicity, the original polynomial coefficients b'^ ( IIV.SP 
have been modified to their more convenient form pm through 
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and furthermore the coefficients x^n are given by 

Xmn = l-m-n- . (IV.12) 

Next, consider the auxihary potential A(y) as part of the gauge field ^'^^AQ(r) 
(]III.56|) which becomes a polynomial A2n{y) of degree 2N and whose differ- 
ential equation is then given by (1III.57bl) with the formal solution being 
specified by equation (1III.59I) . i.e. 

Aiy) = ye'^' ■ r %e--'Py' ■ G,^+,{y') . (IV.13) 

Jy'=y y 

Substituting here the polynomial G2N+i{y) from the result fllV.SP yields for 
the desired polynomial v42n(?/) the following form, simply by straight-forward 
calculation: 

A2n{y)= Yl b'Mm + n + l)\drnniy) ■ (IV.14) 

n,m=0 
(n+m>l) 

Here, the coefficient functions amniy) are given by 

Thus, one can easily convince oneself of the fact that the polynomial A2n{y) 
is of the order A2n{y) = o{y) at the origin {y ^ 0). This verifies explicitly 
the conclusion flIII.60p drawn from the general form of A{y) flIII.591) . And 
furthermore with the polynomial A2m{y) being known, it becomes a simple 
exercise to check also the validity of the second integral identity flIII.66p in 
its polynomial form: 

dye-^^yA{y)Q2n{y) = -mf dy A2n{y)e-^^y . (IV.16) 
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The explicit numerical realization of this identity can be deduced from the 
first identitydlVaOj). 

Finally, the quartic identity flIIL67|) is to be specified down to the poly- 



nomial ansatz which then yields by straight-forward calculation 

dy I '-"''^^--^^y^ ] =- I dyAMy)Q2dy)e-''' 







^ >^ (m + n + l)!(/ + g + l)! 



(IV.17) 



n,m=0 l,q=0 
(n+m>l) (l+q>l} 

where the coefficients Tmn,iq are given by 

Obviously, this is the quartic counterpart of the quadratic identities ( llV.lOp 
and (fIVlHIl . 

But now that all constituents of the field energy e!^ (]III.53p . or its mass 
equivalent M^^'^c^ (IIII.61D . resp., are explicitly known for the polynomial 
ansatz, the final result appears in the following form: 

o2 



= -—/?(! + V^ipo,Pi, . . .Pn)) (IV.19) 
Ob 



with the electric potential function Vn being given by 

N 

Vfiipo, Pi, . . . Pn) = 2 ■ ^ PmPn ■ (m + n)! x^„ 



n,m=0 
(n+m>l) 

A A (m + 72 + l)!(/ + g + l)! 

+ 2^ 2^ PmPnPlPq ' 2™+"+'+'?+2 '^rnn,lq 

n,m=0 l,q=0 
(n+m>l) (l+q>l) 

(IV.20) 
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Fortunately, the potential function emerges here as an ordinary (but 
multi-dimensional) polynomial when written in terms of the modified ansatz 
parameters pm. fllV.llI) . This completes and verifies the identity of the field 
energy E^^ and its mass equivalent M^^^c^, which is one of the two physi- 
cal constituents of the energy functional eI^'^ flIII.22p . And furthermore the 
kinematical constraints are automatically obeyed, too, since both the nor- 
malization condition (]III.30|) for the trial function R{r) is satisfied in the 
form (lIII.Sip and the Poisson constraint (lIII.3ip is satisfied in form of the 
integral identities flIII.62p - (IIII.67p ! Observe, however that the lowest-order 
coefficient feg is to be considered as a function of the other ansatz parame- 
ters (3 and h'^{k = 1 . . .N), see equation fllV.Sp which reads in terms of the 
modified ansatz parameters Pm (lIV-lip : 

N N 
pI+ ^ pmPn- {m + n+ 1)\ = PmPnjm + u + 1)\ = 1 . (IV.21) 

n,m=0 m,n=0 
(n+m>l) 

Evidently, this constraint defines a certain compact subspace of the configu- 
ration space parameterized by the set of ansatz parameters {P;po,pi,...pfq}. 

In order to make complete the non-relativistic energy functional 4°^ 
( IIIL22p . it becomes necessary to specify also the kinetic energy E'kin of both 
particles in terms of the polynomial ansatz. The spherically symmetric ap- 
proximation of i^kin has already been determined in equation flIII.28P : and if 
the polynomial ansatz (jIV.ip for the wave function R{r) (1III.35P is substi- 
tuted therein, one arrives at the following form of the kinetic energy: 

2^kin = —/?'(!+ Tn(po, Pi,... Pn)) , (IV.22) 
Ob 



55 



with the kinetic function Tfq being given by 

N 

^n(Po,Pi, • • -Pn) = -4 • ^ PmPn-m'^{m + n-l)\ (IV.23) 

m,n=0 
(m+n)>l 

Thus the value of the energy functional e!^^ upon the set of hydrogen-like 
wave functions is explicitly known so that one can go to look now for the 
stationary points of the corresponding function E!^\j3;pk] in any subspace of 
order N. 

Observe here for this purpose that the non-relativistic energy functional 
(IIII.22P in its spherically symmetric approximation (1III.32P emerges in the 
following form: 

4°){N} = — (/?2(i + Tn) -/?(! + Vn)) . (IV.24) 
Ob 

This is a very pleasant result because it admits to eliminate the exponen- 
tial ansatz parameter (3 for the extremalization procedure. Recall that this 
procedure consists in nothing else than putting to zero all the partial deriva- 
tives of the energy function 4°"* (HV-24|) with regard of the normalization 
constraint flIV.2ip . Thus the extremalization with respect to (3 yields 

^ = — {2/3(1 + Tn) - (1 + Vn)) = (IV.25) 
ap ttB 

which fixes the parameter j3 in terms of the other parameters po,Pi, • • - Pn 
through 

Substituting this result for P back into the energy function (]IV.24p recasts 
the latter to the form 

4){N} = ■ ii±M . (IV.27) 

^ ^ ^ 4aB 1 + Tn ^ ^ 
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From this result one concludes that the total energy £^r^ {N} is always 
negative because the kinetic energy flIV.22p is always positive. Therefore 
one expects that the RST principle of minimal energy will yield for any N 
the highest energy level which is admitted by the applied variational ansatz. 
Whether an analogous conclusion does hold also with respect to a lowest- 
possible energy level depends upon the existence of a lower bound of the 
energy functional, but this is suggested by the subsequent numerical calcu- 
lations. Observe also that the exponential parameter (3 does not enter the 
normalization constraint flIV.2ip : and therefore the extremalization proce- 
dure is reduced to merely an N-dimensional problem over just that compact 
submanifold which is specified by the normalization constraint (]IV.2ip . 

1. Zero- Order Approximation (N = 0) 

The zero-order configuration (N = 0) naturally represents the lowest- 
order approximation to the positronium groundstate. Since this has already 
been discussed in two preceding papers |3l [7j, it may be sufficient here to 
briefly exemplify the working of the principle of minimal energy in the sim- 
plest case. The point with this zero-order approximation is that it sets the 
overall geometric background for the excited states and yields the exact 
groundstate energy of the conventional theory (11.3^ . Therefore it is espe- 
cially interesting to see whether and to what extent this lowest-order result 
becomes modified subsequently by the higher-order approximations (N > 1). 

The polynomial Pn{y) fllV.ip of zero order is 

Poiy) = bo (= const.) (IV.28) 
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and thus the corresponding trial wave function R{r) (1III.35P looks as follows 



R(r) ^ W^(r) = J^boexpi-f3y) = J ^ b', expi- (3y) 



The normalization condition fllV.Sp then identifies both ansatz parameters 6q 
and P, i.e. 

b', = 2P^Po = l, (IV.30) 
so that the normalization constraint ( 1III.30I) is actually obeyed: 



oo 2 



drr (^^""^Rir)^ = (2(3)^ j dy ye-^^^^^ = I . (IV.31) 
Furthermore, the electric field strength WF(r) flIII.47p becomes 



[plF(r) ^ W^(r) = ■ - ■ F{y) (IV.32) 

Ob ^ y^ 

with the function being given by equation (1111.480 . The reason for this 
simple result is that the field strength constituent G{y) (1111.520 vanishes since 
the function Q2N{y) (1IV.2P is also zero for N = 0: 



Q2Niy) ^Qo = K - W = . (IV.33) 

But clearly, if the electric field strength adopts such a simple form, the corre- 
sponding potential [pUo(r) (IIII.56P must be expected to be of a comparably 
simple shape. However, one obviously has 

[pUo(r) ^°Uo(r) = ^ . M (IV.34) 

flB y 

with the function A{y) being given by equation (IIII.SSp . The reason for 
the simplicity of the groundstate result is again that the constituent A{y) 
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fllll.59p of the electric potential ^P%{r) fllll.56p is zero for the groundstate 
(^2N(y) =^ ^o{y), cf. (IIV.13P ). Observe here that both the groundstate 
potential ^°Uo(r) flIV.34D and its field strength ^°>F(r) flIV.32p are finite at 
the origin {y = 0) 



WAo|_ =2/?-^ (IV.35a) 



1 2^=0 aB 



which verifies the more general conclusion for WF(0) flIII.46p and [pUo(O) 
( ]III.60p . For a sketch of this "Coulomb-like" type of potentials, see ref. [7]. 
But clearly for the groundstate situation with vanishing functions Qo{y) 
(]IV.2p and Gi{y) (lIV.Sp . the quadratic integral identity (IIV.IOP is satisfied 
trivially, as well as the identity (1IV.16P and the quartic identity (lIV.lTp . 

Naturally for such a simple field configuration, the corresponding total 
energy ^^^^{0} flIV.27D must be expected to adopt a very simple form, too: 

4>} = -^.(l±]^ = -^. (,V,36) 

4aB 1 + To 4aB 



since both the potential function (IIV.20P and the kinetic function T^ 
f lIV.23p become zero for N = 



\/o = To = . (IV.37) 



The result ( 1IV.36p . however, is rather amazing because it just coincides with 



the conventional groundstate prediction (1I.3P for the principal quantum num- 
ber rip = 0. Observe here, that the present RST result (1IV.36P is an approx- 
imation due to the applied variational procedure, whereas the conventional 
result fll.3p is usually considered to be exact within the framework of the con- 
ventional non-relativistic quantum mechanics] (See appendix). Therefore 
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it remains to be discussed whether or not the true (but hitherto unknown) 
non-relativistic RST prediction does deviate from the (necessarily relativis- 
tic) experimental data [8J to a smaller or larger extent in comparison to the 
deviations of the conventional (non-relativistic) predictions, see below. 

2. First- Order Spectrum (N = 1) 

Obviously, the general formalism in Sect. Ill singles out the zero-order 
approximation as an especially simple subcase. The reason for this is that the 
parameter (3 is of the exponential type whereas all the other parameters bn 
(IIV.II) are of the polynomial type. But on principle, the variational formal- 
ism works in an analogous manner also for the higher-order approximations 
(N > 1). Therefore it may be sufficient here to collect, for the discussion of 
the first-order spectrum (N = 1), all the prerequisites in a shortened proce- 
dure: 

For N = 1, the kinetic function Tpj flIV.23P is easily seen to adopt the 
following simple shape 

Ti(po,Pi) = -4pi(po + Pi) , (IV.38) 

whereas the potential function (IIV.20p must necessarily look somewhat 
more complicated: 

21 

Vi{Po,Pi) = -Pi(2po + 5pi - 2poPi - 9poPi - yP?) • (IV.39) 



With these results, the first-order energy functional E}^'{1} I WIT} appears 



as 



4°){1} = --^.^1±^ (IV.40) 
^ ^ ^ 4aB 1 + Ti ^ ^ 
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and could now be extremalized under the normalization constraint (lIV.2ip . 
The latter is, for the present first-order situation (N = 1), a simple quadratic 
form of two variables {po,Pi), i-e. 

pI + 4poPi + Qpl = l. (IV.41) 

Such a simple extremalization problem over a compact configuration 
space can of course be easily solved, but it is very instructive to look upon 
this procedure also from a somewhat different viewpoint: Since the subspace 
defined by the constraint (]IV.4ip is topologically equivalent to a circle S^, 
one introduces the angular variable a (0 < a < 27r) through 

sin a . 
pi = (IV.42a) 

sin a (IV. 42b) 



Po = cos a 



so that the normalization constraint ( 1IV.41I) is automatically satisfied. Fur- 
thermore, the kinetic function Ti(po,Pi) flIV.381) and the potential func- 
tion Vi{po,pi) flIV.39p become angular functions Ti{a) and Vi(a) through 
the parameterization flIV.42ap - flIV.42bl) . i.e. 

Ti{pq,pi) =^ T'i(a) = 2 sin a ^sina — v^cosaj (IV.43a) 

(1 \/2 7 \ 
- sin a H — — cos a sin^ a — V2 cos a sin^ a I 
2 4 8 y 

(IV.43b) 

which then applies also to the first-order energy flIV.40p 



4aB 1 + 4ai 
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Thus the energy eigenvalue problem (of first oder, N = 1) is reduced to the 
problem of determining the stationary points of the spectral function Si{a) 
without any constraint, see fig. 1. 

Evidently, the first-order spectral function 6*1(0;) flIV.44p has period vr; 
and in this basic interval (0 < a < vr) there exist four stationary points, 
namely two maxima (la, lb) and two minima (lc,ld). The corresponding 

energies 4°^{i} (w:44|) are the following 

= 1 ~ -6, 8029 . . . [eV] {IV Aba) 

rel. maxima: ^ ^ — ^ _ q 086717499 . . . ~ -0, 590 [eV] 

4aB ' L J 



(IV.45b) 



1, 033319474 . . . ~ -7, 030 [eV] 

(IV.46a) 

1,128194657... ~ -7,675 [eV] . 

(IV.46b) 

Here, the first result Eia (lIV.45aP is clearly the most striking one, be- 
cause it agrees exactly with both the zero-order approximation (1IV.36P and 
the conventional groundstate prediction fll.3|) . where the latter is commonly 
considered to be exact within the framework of the generally accepted non- 
relativistic quantum mechanics. However, the occurrence of two minima and 
two maxima in the first-order spectrum demands an explanation: From the 
physical point of view there should occur for N = 1 just one maximum and 



e2 



Eic = 

4a 



rel. minima: ^ 

E-\A = 

4aR 
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one minimum corresponding to the groundstate and the first excited state, 
as is actually the case with the conventional approach (see Appendix). 

But evidently the present approximation procedure generates stationary 
points which are spurious in the sense that they are not present (not even 
approximately) in the next-higher order of approximation. This viewpoint 
receives support from the subsequent consideration of the second-order spec- 
trum (N = 2) which reproduces the energies Eiy, . . . Eid (lIV.45bp - flIV.46bp 
but yields a complete falsification of the maximum E'la (lIV.45ap . Here 
it should be a matter of course that the increasing orders of approxima- 
tion (N N + 1 N + 2...) must be required to produce smaller and 
smaller corrections to the lower-order predictions. This natural requirement 
is satisfied for the relative minima Ei^ and i?id f lIV.46aP - flIV.46bp which may 
be conceived as first-order successors of the single relative minimum f]IV.36P 
of zero order; and thus the numerical coincidence of the first-order maxi- 
mum E'la flIV.45ap with the zero-order minimum f lIV.36p must appear as an 
incident. As a consequence for the present approximation procedure, one 
is naturally led to a selection criterion, namely in the sense that for any 
approximation order N those stationary points have to be rejected as arte- 
facts which are not (even approximately) reproduced in the next order N+ 1. 
For N ^ oo, the stationary points of infinite order would then represent the 
true RST energy spectrum of non-relativistic positronium (in the spherically 
symmetric approximation). This situation is in contrast to the conventional 
treatment (see Appendix), where in any order N there arise just N+1 station- 
ary points which do reappear in all higher orders (N^N+1 ^N + 2...) 
and therefore do represent exact solutions of the conventional Schrodinger 
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equation (II. Sp . 
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-5j(a) -0,6 




Fig. 1: First-Order Spectrum (N = 1) 

Along the "polynomial valley" (lIV.4ip . which is topologically equivalent 
to the circle S*^ (0 < a < 27r), there occur four different stationary points 
(IIV.45a|) - flIV.46b|) of the energy functional e!^\i} flIV.44p . namely two min- 
ima and two maxima. The maximum (la) must be rejected according to 
the selection criterion because this does not survive (as stationary point) the 
transition to the next order (N = 2) of the spectrum. The two minima (Ic) 
and (Id) represent first-order approximations of the true groundstate which 
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attainable only in an asymptotic sense 
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3. Second- Order Spectrum (N = 2) 

For the second order (N = 2) the general energy functional (1IV.27P be- 
comes specialized to 

with the general kinetic function flIV.231) emerging here as 

T2{po,Pi,P2) = -4pi (po + Pi + IOP2) - IGP2 (Po + 6^2) , (IV.48) 
and analogously the potential function ( ]IV.20p is found as 



^2(po,Pi,P2) = PoPi (-2 + 9 - pi + 2 ■ popi + 18 ■P0P2 + l2Q-pip2 + 627 ■ pj) 

+ P0P2 (-10 + 42 ■ P0P2 + 1095 ■ pI) 

f 2 3057 10845 ^\ 

+ P1P2 i -51 + 201 ■ pI + ^ ■ P1P2 + ■ pIJ 

2 21 4 2 15165 4 

-5-p? + y ■p^-147-p^ + __.p4 ^ 

(IV.49) 

Furthermore, the general normalization condition ( 1IV.21I) adopts now the 
special shape 

pI + 4poPi + 6pi + 120p^ + 12poP2 + 48pip2 = 1 (IV. 50) 

which defines some compact submanifold topologically equivalent to a two- 
sphere S"^; and therefore the following angular parameterization is feasible 



67 



Po = cosai — v^sinai cosa2 + v^sinai sin 02 (IV.51a) 

Pi = sin «! cos ^2 — V 3 sin a\ sin 0^2 (IV. 5 lb) 
a/3 

P2 = — sin «! sin 02 (IV. 51c) 
6 

(0 < ai < vr; < ^2 < 27r) . 

This converts again the kinetic function T2{po,pi,p2) and the potential func- 
tion V2{po,pi,p2) to the corresponding functions over the 2-sphere S"^, i.e. 



72(^0,^1,^2)^72(01,02) (IV.52a) 
V2{po,Pi,P2) ^ V2iai,a2) . (IV.52b) 

The stationary points of the spectral function 5*2 (ai, 0:2) 

e / N • (1 + ^2(ai,a2))^ 

S2{ai,a2) =F , .rpf r (IV.53) 

1 + T2(ai,a2) 

will then determine the energy again without any constraint: 

4°){2} E2{aua2) = --^ ■ 52(ai,«2) ■ (IV.54) 

4aB 

Clearly, this is now the second-order generalization of the first-order result 

The spectral function S'2(ai,a2) displays the occurrence of eight mirror 
pairs of stationary points. This mirror arrangement is a consequence of the 
invariance of the spectral function 5*2(0:1, 02) with respect to the inversion: 

{Po,Pi,P2} {-Po, -Pu -P2} , (IV.55) 
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or expressed in terms of the spherical parameterization {ai, 02) 



«! ^ TT — ai (IV.56a) 
^2 ^ TT + tta • (IV.56b) 

Of course, this invariance of the energy functional e!^^ with respect to the 
reflection fllV.SSp is nothing else than the fact that the wave function can 
be replaced by its negative (\E' => — ^) without changing the physics: both 
the mass eigenvalue equations flII.45aP - flII.45bp and the Poisson equations 
f lll.43ap - flll.43dp are invariant with respect to the inversion flIV.55p . 



For the second order of approximation (N = 2), our numerical program 
finds nine stationary points (apart from reflection), where three of them do 
occur already in the first order (N = 1): these are the configurations (lb) 
f lIV.45bp . (Ic) flIV.46ap . and (Id) flIV.46bp . This is a similar phenomenon as 



did occur for proceeding from the zero-order approximation f]IV.36P to the 
first-order level i?ia f lIV.45ap . However the latter configuration (la) unfortu- 
nately is the only first-order case which is not validated by the second-order 
approach (N = 2), despite the fact that it exactly agrees with the conven- 
tional result! Instead, the six newly emerging configurations for N = 2 are 
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the following: 

E2e = —. 0, 0259932 . . . ~ -0, 177 . . . [eV] {IV. 57a) 

E2f = —. 0, 0925612 . . . ~ -0, 629 . . . [eV] (IV.57b) 

E2g = ■ 0, 1130532 . . . ~ -0, 769 . . . [eV] (IV.57c) 

= 1, 0437688 . . . ~ -7, 100 . . . [eV] (IV.57d) 

4aB 

E2i = —. 1, 1286357 . . . ~ -7, 677 . . . [eV] (IV.57e) 

4aB 

Eaj = 1, 1290162 . . . ~ -7, 680 . . . [eV] (lV.57f) 

4aB 

where the notation is chosen in such a way that -Eib = E2h, Ei^ = E2c, -Eid = 
i?2d- The comparison of this second-order spectrum to the preceding first- 
order spectrum (lIV.45al) - (lIV.46bl) demonstrates that the energy predictions 
are still strongly changing when the approximation order N varies from N = 1 
to N = 2. Thus one concludes that the order N has to be increased consider- 
ably in order to get "stable" results (i.e. with negligible dependence upon N). 
In contrast to this, the conventional predictions are absolutely stable with 
respect to any chosen order N (see Appendix). 

The situation is more favorable with the minimal energies of both cases 
N = 1 and N = 2. It is rather obvious that all five energy values 
{Eic, -Eld, E2h, E2i, E2}) are relatively close to the "exact" value (lIV.36p . This 
suggests the conclusion that there exists a true minimum at (roughly) —7, 7 [eV] 
so that the present N = 1 and N = 2 minima yield a more or less good approx- 
imation to this lowest-energy state. Since the conventional and experimental 
values are in the vicinity of (roughly) —6,803 [eV^], there is a discrepancy 
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of 0, 9 [eV] with respect to the present spherically-symmetric RST predic- 
tions. But this discrepancy is surely not due to the relativistic effects which 
are of order as ^ ~ 10^^ . . . 10~^. Thus it remains to be clarified whether per- 
haps the present discrepancy is a consequence of the spherically symmetric 
approximation (Sect. 111.3) or whether RST is intrinsically inaccurate. But 
here it is interesting to note that a first estimate of the anisotropic effect [3l [7] 
yielded a correction contribution of (roughly) 1 [eV] which is potentially just 
in the right order of magnitude in order to let the "exact" non-relativistic 
RST prediction for the positronium groundstate (practically) coincide with 
the conventional prediction of —6,8 [eV]. 

4- Third- Order Approximation (N = 3) 

On principle, the order of the spectrum could be increased arbitrarily 
provided the computer capacity is large enough. Being satisfied with a certain 
lack of accuracy, our computer program manages the case with N = 3 and 
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finds the following energy levels: 



.2 



= 0, 0112637 . . . ~ -0, 0766 , 

4aB 



^2 



Esb = 0, 0259935 . . . ~ -0, 1768 , 

4aB 



o2 



= 0, 0283261 . . . ~ -0, 1926 , 

4aB 



^2 



= 0, 0382245 . . . ~ -0, 2600 , 

4aB 



,2 



Ese = 0, 0925615 . . . ~ -0, 6294 , 

4aB 

= 0, 1130522 . . . ~ -0, 7690 , 

4aB 

Es„ = 0, 1159734 . . . ~ -0, 7889 , 

4aB 

e2 

= 1, 0437690 . . . ~ -7, 1006 , 

4aB 

^3; = 1, 0484095 . . . ~ -7, 1322 , 

4aB 

= 1, 1286357 . . . ~ -7, 6779 , 



4aB 

„2 



eV] 


(IV.58a 


eV] 


(IV. 58b 


eV] 


(IV.58C 


eV] 


(IV.58d 


eV] 


(IV.58e 


eV] 


(IV.58f 


eV] 


(IV.58g 


eV] 


(IV.5811 


eV] 


(IV.58i 


eV] 


(IV.58j 


eV] 


(IV.58k 


eV] 


(IV.581 



= 1, 1288228 . . . ~ -7, 6792 . 

4aB 

e2 

= 1, 1290162 . . . ~ -7, 6805 . 

4aB 

The general conclusion from such an arrangement is that all the energy levels 
newly emerging in the order N will be present also in the order N + 1, while 
all the levels of the preceding order N — 1 have disappeared (check this also 
for the passage from N = 1 flIV.45aD - flIV.46bD to N = 2 flIV.57aD - flIV.57ep l 
For the present orders N = 2 and N = 3 the common energy levels are the 
following: 
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E2e ^ ^3b ^ 0, 0259935 . . . ~ -0, 176 . . . [eV\ (IV.59a) 

4aB 



^2 



£^2f ^ £^3e - -1 0, 092561 . . . ~ -0, 629 . . . [eV] (IV.59b) 

4aB 



o2 



^2g ^ ^3f ^ --; 0, 11305 . . . ~ -0, 769 . . . [eV] (IV.59c) 

4aB 



o2 



^2h <^ ^3h ^ 1, 04376 . . . ~ -7, 100 . . . [eF] (IV.59d) 

4aB 



■,2 



£;2i <^ £;3j ~ 1, 1286357 . . . ~ -7, 677 . . . [eV\ (IV.59e) 

4aB 



,2 



Eaj ^ £^31 ^ 1, 1290162 . . . ~ -7, 680 . . . [eV] (IV.59f) 

4aB 

For the relative arrangement of these thh^d-order (N = 3) RST results in 
comparison to their conventional counterparts see fig. 2. This arrangement 
(logarithmic scale) seems to support the hypothesis that the RST predictions 
will come even closer to their conventional counterparts when the anisotropy 
of the RST interaction potential is taken into account. It appears also some- 
what amazing that the conventional levels for rip = 1 and rip = 3 should have 
no RST counterpart (at least up to the present order N = 3). 
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Fig. 2: RST Energy Levels of Third Order (N = 3) 

The RST levels (llV.58aP - fllV.581l) are relatively close (but not identical) to 
the conventional levels (0 < np < 5) as must be expected for the spherically 
symmetric approximation. The highest RST level E'sa (]lV.58a|) is inclusive 
because the order N must be increased considerably (N >> 3) in order to 
describe such high energy levels. But the conventional groundstate (rip = 0) 
is relatively good approximated already in the present order N = 3. 
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A Conventional Treatment 



The conventional treatment of the internal Coulomb force problem (11.50 
yields some further insight into the peculiarities of the variational approx- 
imation method of Sect. IV. Such a conventional approach relies upon the 
energy functional -Econv 



— / d'r 

2m 



V0c(r) 



1 2 



(A.l) 



whose variational equation is just the conventional Schrodinger equation 
( 11.51) . provided the constraint of wave function normalization is respected 

rfV0c(f)2 = l. (A.2) 

r) into 

(A.3) 



1 . 

Here it is convenient to recast the conventional wave function 
dimensionless form (^{y): 



1 



with the "reduced" Bohr radius Ob 



ah 



2aT 



(A.4) 



being due to the reduced mass m = M/2; and then the conventional func- 
tional (lA.ip reappears in dimensionless form as 



E = 

-"-^conv 

ah 



2 / dyy^ 



\ dy J 



dyy ■ (Piyf 



(A.5) 



where the original normalization condition (IA.2I) reads now in dimensionless 
form 

/•oo 

dyy^<P{yY = l. (A.6) 
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Obviously, the present energy functional (lA.Sp with (1A.6P is the conven- 
tional analogue of the non-relativistic RST functional e!^^ fllll.321) . But here 
it is important to realize the crucial difference of both approaches: whereas 
the conventional theory describes the mutual interaction by means of the 
rigid Coulomb potential with no intrinsic dynamical freedom, RST adopts 
the interaction potential ^^^Aq as a proper dynamical variable obeying its own 
field equations, even in the non-relativistic approximation (see the Poisson 
equation fllll.341) ). 

But despite this structural difference of both approaches, the set of trial 
functions for extremalizing either functional fllll.32|) or flA.5l) may well be 
taken to be the same, i.e. one may try the same functional form for the 
present conventional treatment: 

My) = V2^Pn{y)e-^y =^(j2 bnyA ^"^^ , (A.7) 

\n=0 J 

as was done for the preceding RST treatment, see (1111. 35|) with (lIV.ip . Defin- 



ing also the modified polynomial coefficients Pm quite analogously to the RST 
case (HVTT]) 

^ (2/5)™+i ' ^^'^"^ 
the conventional normalization condition flA.6l) then transcribes to these co- 
efficients as 

N 

J2 Pn^Pni^ + n + 2)\ = 1 . (A.9) 

m,n=0 

Clearly, this defines again a compact submanifold of the configuration space 
spanned by the hydrogen-like trial functions flA.7l) . 

It should be obvious that the RST and the conventional approach share 
many common features; and this becomes even more evident by considering 
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now the action of the conventional energy functional (lA.Sp upon the set of 
the hydrogen- like wave functions. First, the kinetic part (£^kin, say) of the 
conventional energy -Econv (I A. 51) looks quite similar to the RST case flIV.22p . 
namely 

4in{N} = ^ fdyy'f^^y = ^P'{l + %{po,Pu...Pn)) (A.IO) 
2ab J \ dy J 2a^ 

where the conventional analogue of the RST kinetic function ( 1IV.23P 

is given by 

N 

%i{pQ,pi,. . .pm) = p^Pn{m + n)\ ■ {w? + + n + m} . (A.ll) 

m,n — 

(m+n)>l 

But naturally, the conventional analogue (Vn) of the RST potential func- 
tion V^ flIV.20p will appear to be considerably simpler because, in the conven- 
tional approach, both particles do interact instantaneously via the standard 
Coulomb potential and this results in a purely quadratic potential func- 
tion Vn- More concretely, the potential energy contribution E^^^ to i^^conv 
flA.SP is given by 



If 2 

^pot{N} = -— / rfyy ■ 01,(1/) = -—/3(1 + Vn(po,Pi,...Pn)) , (A.12) 

where the conventional function Vn is now found to be merely a quadratic 
form of the polynomial coefficients Pm'- 

N 

Vn(Po,Pi, • • -Pn) = - ^ PmPn{m + n + 1)\ ■ {m + n) . (A.13) 

m,n— 

(m+n)>l 

The comparison of this result to its RST analogue Vn flIV.201) says that 
the dynamical character of the RST interaction is encoded in form of a quar- 
tic contribution (~ PmPnPiPq) to the potential function V^. Naturally, this 
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structural difference of both approaches will have its consequences for the 
non-relativistic predictions of the positronium level system. The conventional 
description may now be based upon the total energy function i^convlN}, i.e. 

'1 



^conv{N} = 4in{N} + fpot{N} = — 

Ob 



This is the conventional analogue of the RST functional £'.^°^{N} ( IIV.24P : 



. (A. 14) 



but unfortunately it is presently possible to analyze the value of this RST 
functional (1IV.24P only in the spherically symmetric approximation. In con- 
trast to this, the extremalization of the conventional functional i?conv{N} 
flA.14p is easily feasible and corresponds to the ordinary Ritz variational 
principle [31 E] • 

1. Zero- Order Approximation (N = 0) 

The lowest-order polynomial Po{y) (1A.7P just consists of the lowest-order 



coefficient bo (= 2/3po) which, by means of the normalization condition (IA.9p . 
is related to the variational parameter j3 through 

bo = 2Ppo = V2P^po = -^ . (A.15) 



Consequently, the normalized trial function flA.7p of zero order becomes 



My) = ^P'^e-fy , (A.16) 

while both the kinetic and potential functions do vanish trivially {To = Vq = 
0). Thus the corresponding total energy i?conv{0} (1A.14P adopts a very simple 
form: 



i?conv{0} = — (A.17) 



7J 



whose minimal value is given by 

dp '/^-^ 2ab 4aB 

This just coincides with both the conventional Schrodinger result (11.30 and 

the RST prediction e5?^{0} of zero order fllV.36p . 



But clearly, one cannot be satisfied with such a zero-order success; and 
this forces one to consider the next higher order. 

2. First-Order Spectrum (N = 1) 

The first-order wave function </)i {y) is deduced from the general expression 
fOTfl) as 

0i(l/) = V^(6o + Me-^\ (A. 19) 



where the polynomial coefficients 60, hi (or po and pi, resp.) must satisfy the 
normalization condition flA.9p . i.e. for N = 1 



2pi + 12poPi + 24pf = 1 . (A.20) 



This is the conventional analogue of the former RST condition flIV.4ip and 



therefore can also be reparametrized by some angle a (0 < a < 27r): 

Pq = ■ cos a — -y^ ■ sin a (A. 21a) 

sin a / . , X 

pi = , (A.21b) 

cf. the RST analogue hereof fllV.42aD - flIV.42bp . 

As a consequence, both the kinetic and potential function Ti{pq, pi) f lA.llD 
and Vi(po,Pi) ( lA.lSp . i.e. 

^i(Po,Pi) = -8pi{po + 2pi) (A.22a) 

Vi(po,Pi) = -4pi(po + 3pi) , (A.22b) 
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become now functions of that angle a: 

4, . o 



- sm a 
3^ 



-\/3 si 



sm a cos a 



Vi(po,Pi) ^ Vi(a) = ^ sin a cos a , 

v3 

and the same does then apply also to the energy functional Ec. 

^2 



^conv{N} ^ E 

c 



.{1} = - 



-/32(1 + Ti(a))-/3(1 + Vi(a)) 



(A.23a) 
(A.23b) 

. (A.24) 



Here, the same procedure for eliminating the trial parameter (3 can be used 
as was done for the RST case flIV.25p . i.e. the equilibrium value for (3 

1 + Vi 



/5 



(A.25) 



l + Ti 

is substituted back to -E'conv{l} in order to yield the conventional energy 
function Eiia) as 



;i + Vi(a)r 



(A.26) 



Aa-Q 1 + 71(0;) ' 4aB 
Of course, this is again the conventional counterpart of the former RST 
result (1IV.44P : but the conventional spectral function Si{a) (]A.26p . as op- 
posed to its RST counterpart 6*1(0;) (fig.l), displays an important differ- 
ence: whereas the RST function 6*1(0) has two maxima (la, lb) and two 
minima (lc,ld) within the basic interval (0 < o < 27r), the conventional 
function Si{a) flA.26P has only one minimum (lA) and one maximum (IB), 
see fig. A.l. The minimum occurs at oia = (mod vr) and therefore the 
polynomial coefficients po,Pi (lA.21a|) -( 1A.21b|) become 



Pi 



lA 



lA 



, 



(A.27a) 
(A.27b) 



80 



which lets reappear the groundstate wave function 0o(l/) (1A.16P of the zero- 
order approximation (N = 0) also here on the first-order level flA.19|) . 



Clearly, with the minimum value Si{a) = 1 of the spectral function Si{a) 

lA 

one is again led back to the groundstate energy flA.lSp of the zero-order 
spectrum. 

But a new element emerges on the first-order level (N = 1) in form of 
the maximum (IB), see fig. Al, which has the following polynomial coeffi- 
cients po, Pi flA.21al) - flA.21bD : 

V2 



Po 
Pi 



IB 



IB 



V2 
4 



(A.28a) 
(A.28b) 



Furthermore, the associated exponential parameter j3 (IA.25I) becomes 



P 



IB 



(A.29) 



so that the wave function (f)i{y) (1A.19P adopts the usual conventional form 
for the standard Coulomb force problem 



>1B 



V2 y _ 

Li - e 

2 4^ 



with Li(|) denoting the first Laguerre polynomial 



2 



(A.30) 



(A.31) 



Correspondingly, the energy £'conv{l} (1A.26P is then found to coincide as 
expected with the standard result (II. 3p for the first excited state (rzp = 1) of 
the ordinary Coulomb force problem: 



1 



IB 



4aR 22 



16aB 



-1,701... [eV] . 



(A.32) 
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Fig. Al: Conventional First-Order Spectrum (N = 1) 

The conventional spectral function Si{a) ( 1A.26I) has one minimum (lA) 
at c^Ij^a = (modTr) and one maximum (IB) at = | (modvr). This 

situation represents the groundstate (-w lA) and first excited state (-w IB) 
of the conventional theory, see equations (lA.lSp and (lA.32p . 
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3. Second-Order Spectrum (N = 2) 

It is both instructive and advantageous to reassure oneself that also 
for N = 2 the variational process must take place over a compact mani- 
fold (i.e. the two-sphere 5*^). In order to realize this clearly one starts with 
the normalization condition (1A.9|) for N = 2, i.e. 



2pl + 12poPi + 24pl + 48poP2 + 240pm + 720^^ = 1 , (A. 33) 

and satisfies this again by passing over to an angular parametrization; namely 

Po = ■ cos ai — ■ sin ai cos 0:2 + sin ai sin 0^2 (A. 34a) 

sinaiCosa2 2 . . 

Pi = -1= ^smaisma2 (A.34bj 

v6 v3 

sinaisina2 ( k oa \ 

V. = (A.34C) 

(^0 < ai < 7r; < as < 27r) . 

Such a reparametrization converts the kinetic and potential functions 72, V2 
( ]A.11|) and (]A.13|) to the corresponding angular functions 

%{po,Pi,P2) = - (SpoPi + 16pi + 48poP2 + 192pip2 + 576pl) '7^(ai, ^2) 

(A.35a) 

Mpo,Pi,P2) = - (4poPi + 12p? + 24pop2 + 144pip2 + 480p^) V2(ai, 02) • 

(A.35b) 

Correspondingly, the conventional second-order energy function ii^conv{2} 
flA.14p can be recast again to the standard angular-dependent form by simply 
eliminating the exponential parameter f3 

P ( \ {1 + V2{ai,a2)f . 

t2{ai,a2) = — , ^ ^ -1 52(01, 02) • (A.36) 

4aB 1 + 72(01, 02) 4aB 
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The stationary points over the two-sphere S'^ can now be determined again 
without any constraint and it should not come as a surprise that the second- 
order spectral function ^2 flA.361) hast just three stationary points on the 
sphere S^: 



^2 


2A 


= S^ 


S2 


2B 


= Si 






1 


S2 


2C 


~ 32 



1 (A.37a) 

(A.37b) 
(A.37c) 



lA 

_ 1 _ 1 

IB - ^ - 4 

_ 1 

~ 9 ■ 



Obviously the first two cases flA.37ap - flA.37b|) yield nothing else than 
the groundstate (]A.18|) and the first excited state (]A.32p which are already 



present in the zero-order spectrum and in the first-order spectrum. Accord- 
ingly, the stationary points coincide with the former ones (1A.15I) . ( lA.27al) - 
(IA.27bl) and f!A.28al) - f!A.28bl) . i.e. for the groundstate (A) 

Po\2A=MiA=Mo = ^^ (A.38a) 
Pi|2A=Pi|iA = 0' (A-38b) 
and similarly for the first excited state (B) 



.V2 
2 

V2 
4 

P2|2B = 0. (A.39c) 



P0I2B = PoliB = (A.39a) 



P1I2B = Pi|iB = =^~r (A.39b) 



Thus the second-order (N = 2) wave functions for the two lowest-energy 
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states are the same as previously cf. (1A.16P and (1A.30P 



(l>2A{y) = My) = Uy) = 26"^ (A.40a) 

V2 



{y) = <PiB{y) = -^L,{^-)e-'^y . (A.40b) 



But a new element emerges in the second-order spectrum in form of the 
second excited state (2C) flA.37cp . This state is due to the configuration 

Po|2c = ±^' Pi|2C = T^' P2|2c = ±^ (A.41a) 

1 /r' 

3 



which yields the wave function 



( . 2 2 2 



>2^{y) = ^[-\ M--^ + _^2)e-3 (A.42) 



with the corresponding energy (]A.36p 



.2 /iN2 



£^conv|2c = ■ J - -0' 76 [e\/] (A.43) 

cf. the conventional spectrum 0I.3P for rip = 2. 

Summarizing, it should appear evident that the conventional spectrum 
of order N yields just the series of N + 1 "exact" positronium levels (11.30 
with < Tip < N, together with their "exact" Schrodinger wave functions. 
Therefore the applied variational method should be applicable also for calcu- 
lating the non-relativistic positronium levels within the framework of RST. 
However a certain complication arises with the latter approach because the 
RST levels become corrected by any higher approximation (N — > N + 1) 
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whereas the conventional higher-order spectrum leaves uncorrected the lower- 
order result. Prom the geometric viewpoint, this means that any station- 
ary point of the lower-dimensional configuration space remains a station- 
ary point after embedding into the next higher-dimensional configuration 
space (N—>N-|-l—>N-|-2— >...). Or conversely, the restriction of 
the conventional spectral function <Sn(q;i, . . . a^) on the sphere jS"^ to the 
subsphere S^~^ does not generate new stationary points on the subsphere. 
Since this nice property is not shared by the RST counterpart ^^(q;! . • • Q;n) 
it seems desirable to look for a more rigorous mathematical characterization 
of those extraordinary functions being due to both the conventional and RST 
type. 
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